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Abstract. We review recent progress on the zero-range process, a model of 
interacting particles which hop between the sites of a lattice with rates that depend 
on the occupancy of the departure site. We discuss several applications which have 
stimulated interest in the model such as shaken granular gases and network dynamics, 
also we discuss how the model may be used as a coarse-grained description of driven 
phase-separating systems. A useful property of the zero-range process is that the steady 
state has a factorised form. We show how this form enables one to analyse in detail 
condensation transitions, wherein a finite fraction of particles accumulate at a single 
site. We review condensation transitions in homogeneous and heterogeneous systems 
and also summarise recent progress in understanding the dynamics of condensation. We 
then turn to several generalisations which also, under certain specified conditions, share 
the property of a factorised steady state. These include several species of particles; 
hop rates which depend on both the departure and the destination sites; continuous 
masses; parallel discrete-time updating; non-conservation of particles and sites. 
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1. Introduction 

The statistical mechanics of nonequilibrium systems is required to understand the 
macroscopic behaviour of processes occurring throughout physics, chemistry, biology and 
even sociology and economics. Nonequilibrium phenomena are encountered whenever 
systems are relaxing towards an equilibrium steady state and also whenever systems 
are driven i.e. maintained away from equilibrium by external forces. Systems of the 
driven kind, which are the main focus of this work, cannot be described in general by 
equilibrium statistical mechanics. Rather, these systems evolve to a nonequilibrium 
steady state. However, though the statistical mechanics of equilibrium steady states 
is well understood, analogous general principles to guide the study of steady states far 
from equilibrium are only just emerging. 

Considerable understanding of how microscopic interactions influence the 
macroscopic properties of nonequilibrium steady states can be gained from the study 
of interacting particle systems |T| E] however. These are systems defined on a lattice 
on which particles hop from site to site, where the precise definition of the stochastic 
particle dynamics is motivated on physical grounds. Nonequilibrium steady states are 
then constructed by driving a current of particles (a conserved quantity) through the 
system. Models of this kind are known as driven diffusive systems 

The ongoing interest in driven diffusive systems has been sustained by the variety 
of the non-trivial behaviour these models exhibit. Even in one dimension, they can 
undergo phase separation and phase transitions. These transitions may be driven by the 
boundary dynamics or by defects for example, or, in translationally invariant systems, 
they may be accompanied by a spontaneously broken symmetry. Furthermore, one 
observes long-range (i.e. power-law) correlations, not only at criticality, but generically in 
the steady state of driven diffusive systems. Notably, these properties of driven diffusive 
systems in one dimension are absent from their one-dimensional classical equilibrium 
counterparts. 

The purpose of this article is to review recent work on a particular driven diffusive 
system: the zero-range process (ZRP). In the ZRP, particles hop from site to site on a 
lattice with a hop rate which depends, most generally, on the site from which it hops and 
the number of particles at the departure site. Despite being simply stated, it displays 
all of the non-trivial properties mentioned above, but with the additional virtue that 
the steady state is given exactly by a factorised form; this simple form of the steady 
state solution offers an opportunity to analyse these properties exactly. 

The ZRP and generalisations we discuss can be employed to investigate fundamental 
aspects of nonequilibrium statistical mechanics, in that one can address issues such as 
the role of conservation laws, the range of interactions, constraints in the dynamics, 
and disorder, all within the framework of an exactly solvable steady state. It has also 
been widely applied as model for nonequilibrium phenomena such as sandpile dynamics 
and the dynamics of avalanches, granular systems, interface growth, polymer dynamics, 
various transport processes, and glasses [I]. 
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The ZRP was previously reviewed in jlj. Since then this model has stimulated 
considerable interest and a number of authors have contributed to a large body of 
results. Here we aim to summarise this progress, with emphasis on developments in the 
'physics literature'. In particular, a deeper understanding of condensation — a transition 
to a phase in which a single site contains a finite fraction of particles in the system - 
has emerged. It has now become apparent that this condensation transition appears 
in a number of unexpected contexts such as wealth condensation in macroeconomies 
[Sj, jamming in traffic [EJIIj, coalescence in granular systems [HUH], gelation in networks 
[TQ1EJ- Further, we discuss the dynamics of the process by which a condensate emerges, 
which have recently become better understood. Finally a number of generalisations of 
the ZRP which retain the property of a factorised steady state have been identified 
and we review these generalisations here. We also discuss some results that have not 
so far appeared in the physics literature such as the relation between condensation in 
the canonical and grand canonical ensembles ( section |4.3.2|) . dynamics which generate 
different ensembles (section 16.4)1 and the Misanthrope process (section lfi.2j) . 

The layout of the review is as follows. In section|2]we define the model and derive the 
steady state. In section El we provide details of some of the nonequilibrium phenomena 
which have recently been modelled using the ZRP. We discuss the homogeneous ZRP 
(i.e. site-independent hop rates) in section H] and show how the model undergoes 
a condensation transition and we compare the analysis of the condensation in the 
canonical and grand canonical ensembles. We then discuss the coarsening process 
which determines the dynamics of the condensation. We turn to condensation in 
the heterogeneous ZRP (i.e. site-dependent hop rates) in section |3] and show how the 
mechanism of condensation and the coarsening dynamics differ from the homogeneous 
ZRP. In section |B] we discuss several generalisations of the ZRP, all of which are still 
characterised by the exact, factorised steady state. These include generalisations to two 
species of conserved particles (i.e. two conservation laws), to a model in which the hop 
rates depend on the number of particles at both the departure and the target sites, to 
continuous mass variables, to allow arbitrary fractions of the mass at a site to hop, to 
a parallel update mechanism, and to relaxing a conservation law by allowing particle 
number and/or site number to fluctuate. We summarise in section 

2. Definition and steady state 

2.1. Definition 

The zero-range process is a model in which many indistinguishable particles occupy 
sites on a lattice. Each lattice site may contain an integer number of particles and these 
particles hop between neighbouring sites with a rate that depends on the number of 
particles at the site of departure. 

In one dimension, the ZRP is defined on a lattice containing L sites, labelled 
I = 1, . . . , L, and we consider periodic boundary conditions (i.e. site L + 1 = site 
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Figure 1. Mapping between the zero-range process and the asymmetric exclusion 
process 



1). The number of particles at site I is ni, an integer greater than or equal to zero. The 
total number of particles in the system is N and the particle density, p, is given by 

We consider totally asymmetric dynamics (although this is generalised in section 12. Hj) 
such that a particle hops from site / to its nearest neighbour site to the right (site I + 1) 
with rate u(n{). 

In general one can consider any lattice in any dimension, including disordered 
lattices, and it turns out that one can still obtain the exact steady state (see section EH]) - 
For simplicity, we will begin by considering the zero-range process in one dimension. 



2.1.1. Mapping to an asymmetric exclusion process A useful property of the zero-range 
process in one dimension is that it can be mapped onto an exclusion process (i.e. a model 
in which lattice sites are either occupied by a single particle or they are vacant). This 
mapping is illustrated in figure [T] The mapping is constructed by identifying a particle 
configuration in the ZRP with a corresponding configuration of particles in an exclusion 
model. (The mapping is unique up to translations of the exclusion process lattice.) To 
do this, one thinks of particles in the ZRP as vacancies in the exclusion process, and 
sites in the ZRP as occupied sites in the exclusion process. Thus, in figure [TJ site 1 in 
the ZRP becomes particle 1 in the exclusion process. The next three vacancies in the 
exclusion process represent the particles at site 2 in the ZRP and then site 2 itself is 
represented by particle 2 in the exclusion process, and so on. In this way one obtains 
an exclusion model on a lattice containing L + N sites and L particles. 

The exclusion process dynamics are inferred from the way in which configurations 
evolve when the corresponding ZRP configurations evolve under the ZRP dynamics: 
the hop rates in the ZRP, which depend on the number of particles at the departure 
site, become hop rates in the exclusion process which depend on the distance to the 
next particle in front. Thus, depending on the form chosen for u(n), there may be a 



6 



long-range interaction between the particles in the exclusion process. 

We remark that this mapping relies on the preservation of the order of particles 
under the exclusion process dynamics therefore it is only really useful in one dimension. 



2.2. Solution of the steady state 

One of the most important properties of the ZRP is that its steady state is given by a 
simple factorised form. This means that the steady state probability P({n{\) of finding 
the system in a configuration {n{\ = n 1 ,n 2 , . . . ,ul is given by a product of (scalar) 
factors f(ni) — one factor for each site of the system — i.e. 

L 

P{{n l })=Z^ N \[f{n l ), (2) 
1=1 

where -Zl.tv is a normalisation which ensures that the sum of probabilities for all 
configurations containing N particles is equal to one, hence 

{n,} i=l \l=l J 

Here, the 5-function has been introduced to guarantee that we only include those 
configurations containing N particles in the sum. Finally, the factors f(ni) are 
determined by the hop rates: 

n 

/(n)=n«(i)" 1 for n>0, /(0) = 1 . (4) 
i=i 

We now turn to the proof of the steady state (J2J) to (HJ). The first step is to write 
the steady state condition that is satisfied by the probabilities P({ni}). This condition 
balances the probability current due to hops into a particular configuration with the 
probability current due to hops out of the same configuration, hence 

L 

= • • , n,_i+l, n,-l, ■ ■ •) - u(n,)P({n,})] 6(n{) , (5) 

1=1 

where the Heaviside function 9{n) is included to emphasise that site I must be occupied 
for there to be associated hops out of and into the configuration {rii}. The next step is 
to substitute the factorised form (J2J) into (jHJ) and look to equate each term in the sum 
separately, hence 

u(m„i + l)/(n,_i+l)/(7i,-l) = u(7i,)/(n,_i)/(n,) , (6) 
after cancelling common factors. This equation in turn implies that 

/(th-i + 1) f(ni) 

w(72/_i + l)— r- — = u(ni)—. -r = constant , (7) 

f{ni-i) f{ni-l) 

for all values of I. The constant can be set equal to unity without loss of generality, 

hence 

fin,) = ^=r • W 



which is readily iterated to yield (jlj) where we have set /(0) = 1, again without loss of 
generality. This completes the proof of the steady state (J2J) to (Jlj). 



2.2.1. Some useful general results From the three equations (j2j) to Q we can, in 
principle, compute all the steady state properties of the ZRP. We remark here that the 
steady state behaviour is determined by the form of f(n). We could, if we desired, 
choose any form we like for f(n) and then infer the hop rates from 

having rearranged (JHJ). The consequence of this is that any model, with integer site 
variables and a conservation law, which has a factorised steady state can be analysed 
within the framework of a corresponding ZRP — one can use the ZRP to provide a 
complete account of all the possible steady state behaviour in such models. 

It is important to note that f(n), the single-site weight, is distinct from p(n), the 
probability that a given site contains n particles, given by 

P(n) = /(n)%^. (10) 

To obtain this equation we fix the occupation of site one to be n then sum (J2J) over 
allowed occupations of the remaining sites subject to the constraint that the remaining 
number of particles is iV — n 

pin) = P(n, n% . . . n L ) 5 ( n; — (N — n) 

n.2,...,nz, \/=2 

Thus the (^-function constraint in © is important in that it induces correlations between 
sites. 

It is also straightforward to obtain an expression of the mean hop rate (u{n)) where 
(■ • ■) denotes an average in the steady state. In the case of asymmetric dynamics (u(n)) 
is just the particle current. In general, 



(u(n)} = — 



' m,...,n£ Z=l \ 1=1 / 



ill) 



Zl,n 

where we have used (JEJ). Thus the mean hop rate is independent of location in the system 
so it is a conserved quantity. This result is not necessarily obvious in the case where 
the dynamics are symmetric for example, in which case (u{n)) remains finite although 
the current vanishes. 

Another useful exact result is a recursion for the partition function Zl^\ 

N 

Z l ,n = ^2f(n)Z L _ 1:N _ n , (12) 

n=0 
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which can be easily obtained by summing (jlUJI over n. This result is usefully employed 
as an algorithm to be iterated on a computer. Thus it is straightforward to compare 
analytic results with numerics for p(n) or (u(n)) for example. 



2.3. Generalisation to an arbitrary lattice 

In this section we generalise the derivation of the steady state given in the previous 
section to the case where the ZRP is defined on an arbitrary lattice. We use the term 
arbitrary lattice in a very general sense to mean any lattice in any dimension which may 
include heterogeneity in connectivity, or in the hop rates. 

The steady state still factorises if the hop rate from site I to k, Uki(ni), has the 
general form 

u k i(ni) = ui(ni)W k i , (13) 

where U/(n;) is a site-dependent function, giving the total rate at which a particle leaves 
site Z if Z is occupied by rii particles, and W k \ is the probability that the particle hops to 
site k. These probabilities define a stochastic matrix for a single particle moving on a 
finite collection of L (in this context, L is the total number of sites in the system) sites, 
therefore 

xy w =i, (14) 

k 

by conservation of probability. 

The matrix W k i defines an arbitrary connectivity of the underlying lattice: if Wki 
is zero then there is no bond from site I to k across which particles may hop. Moreover, 
if Wki 7^ Wik the hopping is not symmetric, thus the matrix Wki encodes the symmetry 
of the hopping dynamics. Finally we can use Wki or ui(n{) to include heterogeneity in 
the hopping dynamics. Thus the definition of the arbitrary lattice implies these three 
properties (connectivity of the lattice, symmetry and heterogeneity). 

The steady state is still given by (J2J) but with the factors f(ni) — > fi(ni) now 
site-dependent. These factors are given by 

ni 



fi(ni) = I) 



si 



for m>0, /(0) = 1, (15) 

^ ,7/1 / 1 

where the si are the steady state weights of a single random walker which moves on a 
lattice with rates Wki- These weights satisfy 

s i = Yl SkW[k ■ ( 16 ) 

k 

The proof of the steady state follows the same steps as the proof given in the 
previous section. The steady state condition on the probabilities P({n z }) given by © 
is now 



^2 Uk{n k + l)W lk P{- ■ ■ , 7ifc+l, • • • ,n { -l, • • -) 

.k^l 
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e{ni) 



(17) 



where we have used (fPfj) in the second term on the rhs. Again, we substitute the steady 
state form (J2J) (but now with f(ni) replaced by fi(ni)) into this equation and look to 
equate each term in the sum separately, hence 

ui{ni)f k (n k )fi(ni) = ^2u k (n k +l)Wi k f k (n k +l)fi(ni-l) , (18) 

k^l 

for rti > 0. Finally, inserting (fTK|) yields the condition Ijlfijl. 

As an example, consider a one-dimensional chain. One can use the W k i to encode 
partial asymmetry in the hop rates. That is, if a particle hops to the right with rate 
ui+ii{n) = pu(n) and to the left with rate ui-n(n) = qu(n), then this corresponds to 
the choice 

W l+ll =p and WU, = g. (19) 

With this choice, the steady state weights, si, satisfy 

(p + q)si = ps/_i + qs l+1 , (20) 

with the solution si = constant for all I. This constant can be taken equal to one without 
loss of generality. Thus the steady state weights (JH) are unmodified by the degree of 
symmetry of the hopping dynamics. Other lattices where Si = constant are hypercubic 
lattices in any dimension and the fully-connected geometry. A consequence of this is 
that if the steady state factorises for a model in one dimension it will factorise in any 
higher dimension. 

As an example where si is not constant we consider symmetric hopping on a body- 
centred cubic lattice. In this case, each corner site is connected to 6 other corner sites 
and 8 body-centred sites. Each body-centred site is connected to 8 corner sites. Now, 
assuming translational invariance, the steady state weights, s c for corner sites and Sb 

for body-centred sites, satisfy 
g 

Sb = -^s c . (21) 
Thus the steady state weights (|T5j) can be written 

7 



fc(ni) 
for corner sites and 



n-l 

n 



ni 

n 



udi) 



(22) 



(23) 



for body-centred sites. Then the steady state probabilities are 

P({m}) = z~ L ] N J] f c {m) II > ( 24 ) 

zee leB 

where C denotes the set of corner sites and B denotes the set of body-centred sites. 

More complicated si may result from disordered hop rates, for example a one- 
dimensional ZRP with Wi + \ i = pi and Wi-i / = qi has been studied in [HJ H2 UH HI] ■ 
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2.4- Mathematical results 

The ZRP was has been of interest to the mathematical physics community for a long 
time. It was first introduced by Spitzer pQ in 1970 and has subsequently received a lot 
of attention in the mathematical literature. Of particular interest is the hydrodynamic 
limit. Although this is not a limit we focus on, in this section we attempt to give a brief 
overview of some of this work. 

Spitzer £Q obtained the invariant measures for the ZRP for general rates u(n) 
for finite L. Subsequently these measures were shown to be invariant even in the 
case of infinite L j!5j . The hydrodynamic limit (i.e. a continuum limit yielding an 
equation for the time evolution of a coarse-grained particle density field) was obtained, 
for the one-dimensional symmetric ZRP in ^H] , and for the one- dimensional asymmetric 
ZRP in fTf\ ITH] . in which case the hydrodynamic equation assumes the form of a 
nonlinear diffusion equation. The hydrodynamic limit in higher dimensions was given 
in [iniEni- For site-dependent rates Ui, without any dependence on occupation number, 
the existence of invariant measures and the hydrodynamic limit have been proven in 

Steady state fluctuation properties of the homogeneous asymmetric ZRP in one 
dimension were analysed by deriving macroscopic evolution equations for a tagged 
particle [22] and it was shown that the motion of such a particle is determined by the 
characteristic lines of the hydrodynamic equation [2H]. The size of the largest cluster, 
in cases where a finite fraction of particles in the system accumulate at a single site, 
was given in [21], where certain conditions on the hop rates were to given to determine 
whether this fraction is less than one or equal to one (in which case, the system excluding 
the cluster site contains a finite number of particles in the limit of infinite system size). 

This summary is far from exhaustive. In particular we have omitted consequences 
of the hydrodynamic limit, such as the existence of central limit theorems and results for 
large deviations. But most of these results derive from the invariant measure in the limit 
of infinite system size (i.e. L — > oo then t — > oo) whereas our approach is to compute 
the steady state on the finite lattice and only then consider the thermodynamic limit 
(i.e. t —>■ oo then L — ► oo). 

3. Applications 

3.1. Shaken granular gases 

A variety of granular systems can be related to zero-range processes, for example models 
of sandpile dynamics [23121] or models of mass transport through a series of traps |27j . 
Another example, which has received a lot of attention more recently, is models of shaken 
granular gases. These models are based on experiments in which a container is divided 
into L equal compartments by walls, where each wall contains a narrow horizontal slit 
at height h. The container is then mounted on a shaker and filled with N particles e.g. 
plastic balls or sand. See figure [2] for an illustration of the experimental setup. When 
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Figure 2. An illustration of the vertically driven container divided into L = 6 
compartments which are connected by holes at height h. The system has periodic 
boundary conditions, so particles which leave the right-hand compartment to the right 
enter the left-hand compartment from the left, and vice versa. 



the system is shaken vertically, the particles hop from one compartment to another. 

The dynamics of this system clearly resemble those of the zero-range process. To 
make the correspondence explicit, one has to find an appropriate form for the ZRP hop 
rates u(n) motivated by the kinetics of particles within a compartment. Two forms for 
u(n) have been considered [HI EB] and they have been compared, in the context of the 
ZRP, in EES . 

The first of these approaches is due to Eggers |HJ- An expression for u(n) is obtained 
by solving the equations of motion found using the kinetic theory of vertically shaken 
granular gases. For two-dimensional circular discs, the result is 

u{n) = u n 2 e- a{n/N)2 , (25) 

where Uq is a constant (which only sets the time scale), n is the number of particles in 
the compartment and a, given by 

-^(l-e)^(0\ (26) 

depends on the system parameters: g, the acceleration due to gravity; r, the radius 
of the discs; e, the coefficient of restitution; A, the amplitude of the driving; /, the 
frequency of the driving; and Q, the width of each compartment. 

For certain values of a the system is found to evolve to a steady state in which a 
single compartment contains most of the particles. In the thermodynamic limit, iV — > oo 
with L fixed, this corresponds to a phase in which the fraction of particles in a single 
compartment is one. For L = 2 this is a second order transition, from a phase in which 
particles are homogeneously distributed amongst the compartments, and is accompanied 
by a spontaneously broken symmetry. The model was originally introduced for L = 2, 
but has subsequently been studied for L = 3 [SO] and arbitrary L J3HIH21, where it was 
shown the transition is first-order for L > 2. The dynamics of the coarsening leading to 
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the condensed phase was studied in 0. Experiments and simulation indicate that the 
system evolves to its steady state in the condensed phase in two stages: firstly, most of 
the particles cluster in a few compartments, then these clusters disappear one by one 
until only a single cluster remains. The mean cluster size is found to grow like [ln(t)] 1 / 2 . 
The appropriate form of the hop rates for a bidisperse granular gas was derived in [33J; 
this model corresponds to a ZRP with two different types of particles, a generalisation 
of the ZRP we discuss in section |U1 Finally, an asymmetric case where a current of 
particles is supported due to holes connecting compartments at different heights has 
been studied |3~1| . In this case the ZRP description does not yield a factorised steady 
state since the hop rates do not fall into the class discussed in section 12.31 but the 
assumption of factorisation provides a mean-field-type approximation. 

The second approach is due to Lipowski and Droz |2Hj, who propose a simpler 
model than that of Eggers which still captures the essence of the phenomena. They 
consider a hop rate of the form 

U{H) = ^ XP (- To + A(l-n/N) ) ' (27) 
where T Q and A are positive constants, which depend on system parameters. The basis 
of Eggers' derivation of (|25|) was that the effective temperature of a granular system 
decreases when the particle density increases. The form (}2"Tj) is the simplest hop rate 
which reproduces this fact, due to the linear n-dependence in the denominator of the 
exponential. The choice ()27|) then describes dynamics under which a particle is selected 
at random and moved to a randomly chosen neighbouring compartment with probability 
exp(-l/[r + A(l-n/JV)]). 

As in Eggers' model, this model undergoes a transition between a homogeneous 
phase and a phase in which most of the particles occupy a single compartment. Again, 
the transition is second-order for L = 2 [2H1ES] and first-order for L > 2 [3l)l l37|. During 
the coarsening stage of the dynamics the mean cluster size is found to grow like hi(t) 

[2E1 HUES [SUEZ!- 

3.2. Networks 

A network is defined as a set of nodes interconnected by links — mathematically a 
network is simply a graph and the links, which may be directed or undirected, are the 
edges of the graph. The interest within the statistical physics community into networks 
has been with regard to the statistics of the connectivity properties. For example, 
the degree of a node is the number of links attached to it and the degree probability 
distribution may exhibit a 'fat tail' or even a power-law asymptotic behaviour. Simple 
models to generate dynamically such networks and statistics have been introduced and 
studied extensively (for recent reviews see |38l 139]). 

3.2.1. Growing and rewiring networks The dynamics of growing networks is introduced 
through stochastic rules for attaching links to nodes. Basically there are two types of 
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dynamical networks: growing networks where links and nodes are added to the networks 
[4"U] and 're-wiring' or 'equilibrium' networks where the number of nodes and links is 
fixed but links re-attach stochastically from one node to another |4*T| l4"2~] 

In growing networks, power-law degree distributions may be generated through 
preferential attachment. The idea is to add a new node to the network at each time 
step and to attach this node to an existing node selected with probability proportional 
to its degree. It has been shown that the linear dependence of attachment rate on degree 
gives rise to a power-law degree distribution [321 • If the attachment rate grows more 
slowly than linearly in the degree of the node one can generate stretched exponential 
degree distributions; if the attachment rate grows faster than linearly in the degree of 
the node, a condensate node is generated which is connected to nearly all other nodes 

muni. 

Generally, condensation (also referred to as gelation in the network literature) 
occurs when a node captures a finite fraction of the total number of links. As well as 
condensation induced by nonlinear preferential attachment one can have condensation 
induced by heterogeneity jllj. Here nodes have an intrinsic fitness and the rate of 
attachment is proportional to this. The two types of condensation are reminiscent of 
the two types of condensation seen in the ZRP to be discussed in sections l4~2l and 151 The 
connection is explicit in the case of rewiring networks which we review in more detail 
in section 16.2.41 In these networks the nodes correspond to sites of the lattice and the 
dynamics of links becomes equivalent to the dynamics of the particles in a generalised 
ZRP on a fully-connected geometry. 

3.2.2. Directed edges and two-species ZRP One can also consider networks where the 
edges are directed. Then each node is characterised by its in-degree and out-degree 
which are the number of edges pointing into the node and the number of edges pointing 
out of the node respectively. Rewiring dynamics consists of a rate for rewiring the 
outgoing end of an edge and a rate for rewiring the ingoing end of an edge |l2] and 
these rates most generally depend on the in- and out-degree of both the departure and 
destination nodes. 

If we think of ingoing edges and the outgoing edges as two species of particles we 
have a generalisation of the two-species ZRP discussed in section 16.11 Again the model 
is defined on a fully-connected geometry and one obtains a condition for factorisation 
similar in form to the factorisation condition for the two-species ZRP • It should be 
possible to analyse condensation transitions within this framework 

3.3. Coarse-grained descriptions 

So far we have presented applications of the ZRP which essentially involve mapping 
various models on to the ZRP at the level of exchange of particles between sites of the 
system. More generally, however, one may think of the sites of the ZRP as representing 
domains of some driven system — this is a most natural picture within the exclusion 
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process interpretation of the ZRP (figure The domains may have some internal 
structure, for example further degrees of freedom, but this all integrated out, and one 
is left with an effective dynamics of exchange of length between domains given by u(n). 



3.3.1. Bus route model An early example of the use of the ZRP as an effective coarse- 
grained description of more complicated microscopic dynamics is the 'bus route model' 
|4"H] . The model is defined on a Id lattice. Each site (bus-stop) is either empty, contains 
a bus (a conserved particle) or contains a passenger (non-conserved quantity). The 
dynamical processes are that passengers arrive at an empty site with rate A; a bus 
moves forward to the next stop with rate 1 if that stop is empty; if the next stop 
contains passengers the bus moves forward with rate (3 and removes the passengers. 
The bus route model can be related to the ZRP by an approximation of a mean-field 
nature in which we integrate out the non-conserved quantity (passengers). The bus- 
route dynamics can be thought of as an exclusion process which is mapped onto ZRP 
dynamics as in figure ^ the buses correspond to the sites of the ZRP and the number of 
bus stops between a bus and the next bus along the route, corresponds to the number 
of particles at the site of the ZRP. The idea is that the hop rate of the buses becomes 
a function of the distance to the next bus ahead. Now, the mean time elapsed since a 
bus-stop was last visited is given by n/v where n is the distance to the next bus ahead 
and v is the steady-state speed. Therefore the mean-field probability that the site next 
to a bus is not occupied by a passenger is exp(— Xn/v). From this probability an effective 
hop rate for a bus into a gap of size n is obtained by averaging the two possible hop 
rates 1, (3: 



This example serves to illustrate how local dynamics may generate an effective ZRP 
dynamics defined by u(n). 

The phenomenon of condensation, wherein a single site of the ZRP contains a finite 
fraction of the density (to be discussed in I4.2j) . would correspond to a finite fraction of 
the bus stops being between two buses i.e. there is a jam of buses in the steady state 
and all the buses arrive at a bus stop at once! However it turns out that since u{n) 
decays exponentially the condition for a strict phase transition in the thermodynamic 
limit is not met, unless we take the passenger arrival rate A — > 0. However on any finite 
system for A sufficiently small, an apparent condensation will be seen. 

In [IS] a dual model to the bus-route model was used to model clogging in pipes. 
Also a model of ant trails introduced in jl^l m ay be mapped exactly onto the bus route 
model. 

3.3.2. Phase separation in one- dimensional systems Recently, using the ZRP as 
the effective description has allowed insight into the long-standing question of phase 
separation in one- dimensional driven systems |4"7| . Within this description phase 
separation is manifested by the emergence of one large domain and this corresponds 




(28) 
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to the phenomenon of condensation in the ZRP. The criterion for phase separation 
corresponds to the criterion for condensation within the ZRP |4"T1 |4"5] . 

The idea is that domains of particles on a one-dimensional lattice can be represented 
by the sites of a ZRP and the dynamics by which domains exchange length can be 
described by u(n). Thus, in this context, u(n) corresponds to a conserved current flowing 
out of a domain of size n and the assumption is that domains are uncorrelated in that 
the current only depends on the length of a domain. The approach reduces a many- 
domain problem to the properties of a single domain; the key step is the identification 
of the effective current u(n). 

As we shall see in section 14.41 the criterion to have condensation is that either the 
hop rates u(n) in the ZRP should decrease to as n increases or else u(n) should decay 
more slowly to an asymptotic value (3 than (3(1 + 2/n). Thus, the criterion for phase 
separation to occur is that the current of particles out of a domain of size n should 
either decay to zero with domain size or decay to a non zero value more slowly than 
(3(1 + 2/n). 

In the development of this approach particular attention has focussed on a class 
of exclusion models with positive particles, negative particles and vacancies |49j . The 
model we consider is defined on a one-dimensional ring of L sites. Each site i is associated 
with a 'spin' variable A site can either be vacant (s« = 0) or occupied by a positive 
(si = +1) or a negative (si = —1) particle. Particles are subject to hard-core repulsion 
and a nearest-neighbour 'ferromagnetic' interaction, defined by the potential 

V = -|J>s m . (29) 

i 

Here < e < 1 is the interaction strength, and the summation runs over all lattice sites. 
The model evolves according to the nearest-neighbour exchange rates 

+ - + — +~ +0-^0+ 0—^-0,(30) 

g(l+AV) 

where AV^ is the difference in the potential V between the initial and final states. One 
can think of this as a generalisation of the Katz-Lebowitz-Spohn driven lattice gas (KLS 
model) to be discussed in section 16.2.21 

The sites of the exclusion process correspond to the vacancies in this model and the 
domains to the consecutive sequences of positive and negative particles as illustrated in 
figure 01 Note that although a domain comprises both positive and negative particles, 
the effective description is in terms of a one-species ZRP. In the case q > 1 (where 
positive particles drift preferentially to the left of negative particles) the steady state 
comprises three pure regions of vacancies, positive particles and negative particles. This 
corresponds to strong condensation where one site of the ZRP holds a fraction tending 
to one of the particles. The current of particles is exponentially small in the system 
size jSU] and u(n) decays exponentially to zero. A more subtle question is whether 
condensation and phase separation can occur for q < 1. 

In the case e = the model reduces to that studied in [2H E2| and an exact solution 
is available via a matrix product ansatz 53J. The mapping of the steady state of the 
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Figure 3. A microscopic configuration of a two-species driven model (bottom) and 
its corresponding configuration in the ZRP (top). 



exclusion process to the steady state of a ZRP becomes exact and allows an exact 
identification of the effective current u(n). In the case of equal densities of positive 
and negative particles the function u(n) is precisely equal to the current flowing in the 
steady state of an exclusion process of length n with open boundary conditions jUj. 
This gives a simple recipe for determining u(n) and provides a numerical test for phase 
separation, even when the open system cannot be solved exactly In the case of 

nonequal densities, the current is given by an exclusion process of length n with periodic 
boundary conditions but with non conservation of the charge of particles jlH|. In either 
case the current decays with domain size as u(n) = w(oo)(l + 3/(2n)) which does not 
satisfy the condensation criterion. However one may have very strong finite size effects 
leading to sharp crossover effects which may give the impression of a phase transition 
on finite systems [53J E3] . 

More generally, when e ^ 0, and the ferromagnetic interaction is switched on, no 
exact solution is available. However using the identification of u{n) suggested by the 
e = case predicts that phase separation (for q < 1) should occur when e is sufficiently 
large |5fi| . This effect is enhanced when the densities of positive and negative particles 
are unequal [IB]- Although the mapping to a ZRP is not expected to be exact when 
e 7^ 0, its predictions have proven accurate so far. 



4. Homogeneous system 

We now return to the ZRP defined in section|21 In this section we study the homogeneous 
system in which hopping rates u(n) and therefore steady-state weights f(n) are site 
independent (see section We will focus on the phenomenon of condensation. As 

we shall see the condition for this to occur depends on the asymptotic behaviour of f(n) 
and therefore u(n). We also compare the analysis of condensation in the canonical and 
grand canonical ensembles. 

As we have seen in the introduction and section El condensation is observed in a 
variety of physical contexts. The factorisation property allows us to analyse exactly 
the condensation mechanism within the ZRP. The first analysis of this kind on a 
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homogeneous system was carried out in j^Zj where the grand canonical partition function 
of a 'balls-in-boxes' model is analysed. It should be noted that condensation also occurs 
in other systems that cannot be solved exactly and a factorised steady state becomes 
an approximation within which to study condensation p)%l I59j. 

Ideally one would wish to demonstrate condensation by analysing the distribution 
of the number of particles p(n) given by (jlOj) . This expression involves Zl,n, defined 
by (j3J), which is the equivalent of the canonical partition function since the number 
of particles is held fixed. However, to understand the condensation phenomenon, it 
is simplest to work within the grand canonical ensemble where the particle number 
fluctuates. 



4-1. Grand canonical ensemble 

We define the grand canonical partition function as 

Z L ( Z ) = Y,Z n ZL,n, (31) 
n=0 

where z is the fugacity which is chosen to fix the density p through 

"'l-^t 1 - < 32) 

Now, using (jSJ), we have 

00 L 

Z L (z) = ** m ]If(™i) = [F(z)] L (33) 

{m i= 0} 1=1 



where 



and 



F(z) = j> m /(™) ( 34 ) 



m=0 



F'(z) , \ 

One views this condition either as defining explicitly p(z), or as defining implicitly z(p). 
The distribution of the number of particles at a given site becomes 

l \ n f( \ Z l-i z n f(n) 
p{n) = z f{n)—— = (36) 

AL f \Z) 

in the grand canonical ensemble. We can then calculate the mean hop rate (u(n)) = 
J2 n p(n)u(n) as 

(u(n)) = z , (37) 

endowing the fugacity with a physical interpretation. We emphasise that here the angle 
brackets denote a steady state average within the grand canonical ensemble. This is to 
be compared with the expression of the mean hop rate in the canonical ensemble, given 
by (HU). 
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/(») ~ ^ZT • (39) 



4-2. Condensation 

It is easy to show, by taking the derivative with respect to z, that the right hand side of 
(jHSJ) is an increasing function of z. Also let the radius of convergence of F(z), defined 
in (jMJ), be z = (3. Thus increasing the value of z corresponds to increasing values of 
the density p(z) until z takes its maximum value (3 which corresponds to the maximal 
value of the density p c 

* = "3$- <*> 

If p c is infinite then for any finite density one can find a solution of ()35j) for z. On the 
other hand, if p c is finite then for p > p c one can no longer satisfy (|35|) and we have 
condensation. 

Let us be specific and consider f(n) with large n asymptotic form 
A 

/3 n n b 

Then z = (3 is clearly the radius of convergence of (f3"""j) an d whether p c is finite or not is 
controlled by F'(f3) which converges for b > 2. Thus for b > 2 we will have condensation. 

To understand what is happening we look at the single-site distribution p(n) fllOl) 
for the number of particles. For n>l this becomes, using the asymptotic form (13*91 . 

^/l) w A 
F(z) n b 

Consider first b < 2 in which case p c = oo. The distribution ()4()j) is a power law 
~ l/n b with an exponential prefactor cutting off the distribution at n ~ l/ln(z//3). 
Increasing values of z imply that the distribution is cut off at larger values of n and 
thus correspond to increasing values of the density. If b < 2 one can obtain any desired 
density p = J2 n n P( n ) by choosing z to cut off the power law at suitable n. 

For b > 2, on the other hand, the power law l/n b has a finite mean. Thus letting the 
cut-off tend to infinity (i.e. z — > f$) only corresponds to a finite density. Therefore for 
p > p c one requires an extra piece to the distribution which represents the condensate. 
The shortfall in the density, p — p c , will be made up by the presence of this piece. 

To confirm that the condensate corresponds to a single site we consider its canonical 
weight: if there is a single condensate site the weight will be Lf(N — N c ) = 
whereas if the condensate is spread over two sites the weight will be L) Y^' n f(N — N c — 
n)f{n) = 0(L 3 ' 2b ), where the prime indicates a sum over terms n = O(N). The factors 
L and (^) come from the number of ways of choosing the condensate site(s). Now the 
condition for condensation, b > 2, implies that 1 — b > 3 — 26, therefore in the limit 
L — > oo configurations with a single condensate are dominant over configurations in 
which the condensate is spread over two (or more) sites. Note that in the condensed 
phase we have used canonical weights, because in this phase the equivalence between 
ensembles breaks down as we discuss in the next subsection. A rigorous proof that the 
condensate corresponds to a single site is given in [BO] . 



V{n) * ■ (40) 
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Figure 4. ln-ln plot of the single-site distribution p(n) v. particle number n. The data 
are obtained by iterating the recursion relation (|12Jl for Zl.tv for L = 1000 and b = 5. 
The circles represent p = 1/4 where the system is in the fluid phase; the diamonds 
represent p = 4 where the system is in the condensed phase. 



4-3. Relation between canonical and grand canonical partition functions 

We now return to the calculation of the canonical partition function and the question 
of equivalence between canonical and grand canonical ensembles. 
We may write using Cauchy's theorem as 

Z L , N = $£L Z L {s) = [F(s)] L , (41) 

where the integral is around a closed contour about the origin in the complex s plane. 
For large N, L (j4~Tj) is dominated by the saddle point of the integral which we denote by 
s = z. Defining 

<p(s) = -p\ns + \n{F(s)] (42) 
the saddle point z is given by <f>'(z) = and this recovers precisely (f55|). Moreover 

Zln * (2,L)^ mz) ^— (43) 

~ (27^)1/2^)11/2 Z N + 1 ( 4 ) 
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Thus when one has a valid saddle point the canonical and grand canonical ensembles 
are equivalent (iTj. However, for a valid saddle point the value of z must be within the 
radius of convergence of F(z) since the integral contour cannot be extended to enclose 
the singularity at s — (3. Thus for p > p c the asymptotics of the canonical partition 
function must take on an alternative form, signalling the phase transition. Only recently 
has a detailed analysis of this form and the condensation mechanism been carried out 
|62j : we give a brief discussion here and refer the reader to jH2]. Also note that the 
equivalence of ensembles has been studied in [60J. There it was shown that for p > p c 
the fluid part of the canonical distribution coincided with the critical grand canonical 
distribution. 



4-3.1. Condensation within canonical ensemble As p / p c the saddle point z 
approaches a branch cut singularity at z — (3. For p > p c there is no valid saddle 
point, but the integral (|41jl is still dominated by the region s ~ (3. Thus in the critical 
region and in the whole of the condensed phase one can expand about s = (3. For f(n) 
of the form (|39|) we let s — (3(1 — u) and expand F(s) for u small as 

*W - «)) = E { -^ Fik) + Bu^ + . . . (45) 

k=0 

where r is the integer part of b and = F(s)\ s=/3 (see e.g. |63j). The second term 
in (J4"5j) is the leading singular part and the coefficient B is determined through 

u r+1 ~ b d r F 

B = lim ^ V- (46) 

«-+o (b - 1) • • • (b - r) du r y 1 

and it can be shown that 

B = AT(1 - b) , (47) 

where T(x) is the usual gamma function. 

In the condensed phase the contour integral in ()41|) can be deformed to run along 
the imaginary axis in the complex u plane. Then the asymptotic behaviour is controlled 
by the nonanalytic term AT(1 — 6)m 6_1 and one obtains 

F{(3) L r°° du 

The asymptotics of this integral can be evaluated by wrapping the contour around the 
negative real axis and one finds that the leading contribution is 

Z l ,n>n c ^ — ^ — (p - p fL h - 1 aU noninte S er b > 2 • ( 49 ) 

Thus we have a distinct supercritical form for the canonical partition function ()49j1 to 
be compared with the grand canonical form ()44|) which holds in the subcritical regime. 
These forms can be used to analyse the piece of p(n) which represents the condensate 
and one finds distinct universal forms for 2 < b < 3 and b > 3 [62 . 



Z l ,n>n c I — exp L 



b -iAT{l-b) 
u(p-p c ) + ... + u b j— + 
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4-3.2. Condensation within grand canonical ensemble In contrast, an analysis of the 
condensate within the grand canonical ensemble requires allowing z > (3, thus one has 
to impose cut-offs on the sums jBlj and modify (p£5j) to 

N 

Z L (z) = [F N (z)} L where F N (z) = f(n)z n . (50) 

n=0 

Here the cut-off N ensures that no site contains more than N, the total number of 
particles in the canonical ensemble, but where the number of particles still fluctuates. 
Then to satisfy (|38|) for p > p c one takes z = (3(1 + rj) with 77 small and chooses rj so 
that pF'{/3) < zF' N (z) = 0(1). The result is 



, = «!+,), ,= (6-2)^ + !^ +0 (I), (51) 

where the coefficient of the 0(1/ N) term has to be chosen to ensure that (|38|) is satisfied. 
This result implies that the mean hopping rate, (u) = z in the grand canonical ensemble, 
overshoots its thermodynamic limit value (3 on a large but finite system [HB] . 

However the condensate is not correctly reproduced by the introduction of a cut-off. 
This is seen by noting that (jHHjl becomes 

p(n) = exp(n?7 — b Inn) (52) 
b \z) 

which has a minimum at n = b/r] and increases to a boundary maximum at the 
maximum allowed value n = N, see figure 03 The piece of p(n) to the right of the 
minimum corresponds to the condensate. However this piece is not centred about the 
value (p — p c )L as the condensate should be. This is because we have to take z > 1 
in order to account for the condensate which implies p(n) must increase at large n 
rather than having a maximum (as is seen in figure HJ. This is in contrast to usual 
Bose condensation, in momentum space in the ideal Bose gas, where the condensate is 
accounted for by taking z / 1 in the appropriate way [61. To correctly describe the 
condensate in the present case one has to work within the canonical ensemble |62j . 

4-4- Condensation condition on hop rates u(n) 

So far our discussion of the condensation mechanism has centred on the single-site weight 
function f(n). In most applications of ZRP, however, one is given the hop rate u(n) (see 
section OJ, therefore we now consider the conditions for condensation implied on u(n). 

For condensation we require that the infinite series zF'(z) = *Y^ = inz n f(ri) 
converges for some value of z. The ratio of successive terms n — 1, n in the series 
is given by 

(n-l)f(n-l) u(n)( 1 _l\ (53) 



znf(n) z \ 

The ratio test [66J tells us that for convergence of the series this ratio should decay more 
slowly than 1 + 1/n, thus for condensation to occur, i.e. convergence at z — f3, u(ri) 
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Figure 5. ln-ln plot of the single-site distribution p(n) given by l|52|) v. particle number 
n. The curve shown is for A/ F(z) — 1, rj — 0.002 and b = 5, and should be compared 
with the diamonds in figure 01 



should decay more slowly than 
u(n) ~ (3 ( 1 + - 



) 



for n^$> 1 . 



(54) 



This corresponds to the asymptotic expression for f{n) (jHSJ) studied in the previous 
subsections. 

If u(n) decays, for example, as 



where < A < 1, one obtains condensation but with the fluid phase now distributed 
according to a stretched exponential distribution p(n) ~ exp(— an 1 ~ x /(l — A)). 

If u(n) increases as n —>■ oo the series has infinite radius of convergence and so 
condensation never occurs. On the other hand if u(n) — > as n — > oo, the series will 
have zero radius of convergence. This implies condensation at any density. Moreover 
the density in the fluid phase will tend to zero, therefore the fraction of particles in the 
condensate will tend to one. 

These conditions on u{n) lead directly to the criteria for phase separation discussed 
in section 13.3.21 When the current of particles out of a domain of size n decays to an 
asymptotic value /?, for large n, more slowly than + 2/n), phase separation occurs 




for n ^> 1 , 



(55) 
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above some critical particle density. This phase separated state consists of a single 
domain containing a finite fraction (less than one) of particles in the system. On the 
other hand, when the current of particles out of a domain of size n decays asymptotically 
to zero, then strong phase separation occurs leading to a state in which a single domain 
contains a fraction of particles equal to one — the remainder of the system contains a 
finite number of particles. 

One can also infer the condensation phenomena that occur in the shaken granular 
gases discussed in section IH.ll In this case, both the hop rates (jZHj) and (j2Tj) decay to 
zero as n — > oo. Thus, in the limit N, L — > oo with p fixed, both models will evolve 
to a steady state in which a single compartment contains a fraction of particles equal 
to one and the remainder of the system contains a finite number of particles. The 
corresponding analysis in the limit iV — > oo with L fixed has been carried out in 



4- 4-1- Explicitly solvable case It turns out that if one takes u(n) = (3(1 + b/n) Vn > 
one can calculate the critical density, given by (JHHj). exactly [001 Ell- In this case one 
has 

fw-rnih-^wnz (56) 

where 

(a) n = a(a + 1) . . . (a + n — 1) , (57) 
is the Pochhammer symbol. Thus 

> + l)(n+l)! 1 ^V /mn (2) B (2) n 



0|5 (b + l) n+1 p(b + l)£f (b + 2) n n\ 



and one finds using the hypergeometric function identity 

y^ (a)n(&)n _ T(c)T(c-a-b) c>a + b m 

^ (c) B n! -r(c-a)r(c-6) C>a + & (58) 

that F(/5) = 6/(6 - 1) and = b/((3(b - 1)(6 - 2)). Thus © yields 

* = ^2 • (59) 

Note however that one does not expect p c to be universal i.e. it generally depends on 
the form of f(n) for small values of n, not just on the asymptotic behaviour. 

4-5. Dynamics of condensation 

In this section we turn to the dynamics of the condensation process. Starting from a 
homogeneous distribution of particles, the condensate emerges as the result of a non- 
trivial coarsening process [IS]. The excess particles accumulate on a number of sites 



24 



and the coarsening is then determined by the exchange of particles between these sites; 
this results in their elimination and hence growth in the mean particle number at such 
sites. The object of this section is to determine the scaling satisfied by this growth. 
We consider the case where the hops rates have the asymptotic form 

u{n) = l + b/n, (60) 

when the system evolves from an initially random distribution of particles at a 
supercritical density p > p c . Two complementary approaches have been taken to analyse 
the dynamics for this case; one is a mean-field approach, and the other is a heuristic 
scaling approach based on random walk arguments. 

4-5.1. Mean- field dynamics The mean-field approach is due to Godreche [HZj and is 
based on earlier studies of the dynamics of urn models [69, 70 . Here we outline this 
approach, and refer the reader to [67J for details. 

In the mean-field system we assume the steady state factorises, and allow particles 
to hop from any site to any other in the system. Thus the probability pt{n) that a site 
contains n particles at time t satisfies a Master equation given by 

_1— = u(n+l)p t (n+l) + (u t )p t (n-l) - u(n)p t (n) - (u t )p t (n) , (61) 

for n > 0, where (ut) = X!m=i M ( m )P*( m ) * s the expectation value of the hop rate at 
time t, and 

^ = u(l) ft (l) - (u t )p t (0) . (62) 

These equations are non-linear, due to the dependence of Ut of Pt{n), however the work 
of jfiTl inUl ED] has shown how one can obtain the late-time dynamics in the condensed 
phase (and also at criticality, although we do not address this point further here). One 
finds that the system decomposes into two distinct parts: one which includes sites that 
have effectively reached a steady state at density p c , and the other which includes sites 
that contain the excess (p — p c )L particles that form the condensate. 

In the condensed phase, in the limit t — > oo, (u t ) — > (u), the steady state 
expectation value of the hop rate. In the steady state, (u) is equal to the fugacity 
z, and in the condensed phase z is given by its maximum allowed value, (3, determined 
by the radius of convergence of the sum (J2U). For the rate (jBllj) [3 = 1 therefore (u) = 1. 
Therefore for late times one sets 

(ut) ~1+Aet, (63) 

where A is an amplitude and is a small time scale, both to be determined later. Next, 
one observes that the terms in (|61j) which reduce the number of particles at a site are 
those which depend on u(n), and those which increase the number of particles at a 
site depend on (u t ). Sites at which n is finite reach an effective steady state and the 
dynamics of the condensation is determined by those sites at which n is large (we will 
refer to sites at which n is extensive as condensate sites). For sufficiently late times 
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these sites are in the scaling regime. Thus one defines the scaling variable x = ne t} and 
the probabilities pt{n) are rewritten in terms of a scaling function g(x) i.e. 

p t {n) ~ e 2 t g(x) . (64) 

Now one substitutes the three equations (|H0| . (jQBl) and into the Master equation 
(jnij, which yields the differential equation 

g"(x) +(^-A + ^j g\x) + (l - A) <?(z) = , (65) 

satisfied by g{x). As a byproduct one also obtains e t ~ ef, hence 

e t ~ t- 1/2 . (66) 

The amplitude A is determined by the properties of g(x) for small and large x and 
is found to be universal in the sense that it only depends on b, therefore the scaling 
function g(x) is also universal in this sense. Properties of the solution of for g(x) 
can be found in [OH EH] • 

The picture that emerges from this mean-field analysis is that the typical number 
of particles at condensate sites scales as t 1 / 2 . Since the total number of particles at such 
sites is (p — p c )L, the number of condensate sites must scale as (p — p c )L -1 / 2 . Hence 
the mean condensate size (m(t)), defined as the total number of particles at condensate 
sites divided by the number of condensate sites at time t, obeys a scaling law 

(m(t)) ~ t 5 , (67) 

defining the exponent 5, where 

5=1/2. (68) 

This result is exact for the ZRP in an infinite number of dimensions, and is also expected 
to hold above an upper critical dimension. Significantly, numerics indicate this is the 
correct result in one dimension for asymmetric hopping dynamics. However, as discussed 
in [HZ] , below the upper critical dimension both dimensionality and bias in the hopping 
dynamics may be expected to play an important role, and numerics indicate that for 
symmetric dynamics the exponent changes to 5 = 1/3. In order to address these issues 
further we turn to an alternative, heuristic scaling approach. 

4-5.2. Random walk argument This approach, given in [00!, * s based on the observation 
(also observed in the mean-field approach) that the dynamics can be divided into distinct 
regimes: (i) a nucleation regime, during which the excess (p — p c )L particles accumulate 
at a finite number of condensate sites such that each condensate site contains a number 
of particles of order L — the remaining sites, which we refer to as bulk sites, have 
converged to the steady state with density p c ; (ii) a coarsening regime, during which 
the condensate sites exchange particles through the bulk - - the bulk is viewed as 
a homogeneous background during this process. The exchange of particles between 
condensate sites leads to the growth of larger condensates at the expense of smaller 
ones. This, in turn, leads to a decrease in the number of condensate sites and an 
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increase in the mean condensate size (m(t)). Thus the coarsening regime persists until 
only a single condensate site remains. The random walk argument yields a prediction of 
the exponent 5 appearing in the scaling law ()67j) which determines the growth of (m(t)), 
the mean condensate size. 

We now outline how the random walk argument is applied to the ZRP in one 
dimension to obtain predictions for the exponent 5 for both asymmetric dynamics and 
symmetric dynamics. During the coarsening regime, the mean rate at which particles 
hop from bulk sites is given by (u) = 1. However at condensate sites particles escape 
with a rate u(n) — (u) = b/n, where n = 0[(p — p c )L\. Therefore a particle escapes from 
a condensate site in a time that scales as O(L). We must now find the time it takes 
such a particle to reach the next condensate site. This time depends on the symmetry 
of the hopping dynamics. 

a) Asymmetric dynamics: In this case, the particle reaches the next condensate site 
to the right with certainty. The particle mobility is given by the mean hop rate 
from occupied sites in the bulk, which is a finite constant given by (u)/(l — p(0)). 
The average distance between two condensate sites scales as O(L) therefore the 
typical time a particle spends in the bulk, having escaped from a condensate site, 
scales as O(L). This time is of the same order in L as the typical escape time from 
condensate sites therefore there are only a finite number of excess particles in the 
bulk at any time. Thus the transit time between two condensate sites does not limit 
the dynamics of the coarsening in any way. The growth of the mean condensate size 
is determined only by the time it takes a condensate site to lose all of its particles. 
This time scales as 0(L x L) therefore the normalised mean condensate size grows 



where the time scale for the coarsening regime, r ~ 0(L 2 ). Therefore the prediction 
5 = 1/2 for asymmetric dynamics is the same as for the mean-field analysis. 

b) Symmetric dynamics: In this case, having left a condensate site, a particle 
performs a symmetric random walk with a diffusion constant given by (u) /(l—p(0)). 
Now, the probability that a particle reaches the next condensate site, which is 
determined by the solution of the 'gamblers' ruin' problem [21], is proportional 
to the inverse separation of condensate sites. Thus this probability is of order 
L, so particles almost certainly return to the condensate site they have just left. 
The typical time it takes a particle to escape from a condensate site is therefore 
0(L 2 ). The transit time for such a particle is given by the first passage time for 
a symmetric random walk which is proportional to the square of the distance i.e. 
0(L 2 ). Therefore, as in the asymmetric case, this time is of the same order in L as 
the typical escape time from condensate sites so the growth of (m(t)) is determined 
only by the time it takes a condensate site to lose all of its particles. This time 
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(70) 



(P~Pc)L \r) 



where the coarsening time scale r ~ 0(L 3 ). Thus 6 = 1/3 for symmetric dynamics. 

The results of simulations are consistent with the predictions of the random walk 
argument in both asymmetric and symmetric cases fof\ If"""] . 

The random walk argument can also be applied to the ZRP in more than one 
dimension [72,. The condensate sites may be viewed as 'trapping' sites amongst 
which the excess particles are exchanged. The key feature of random walks in higher 
dimensions is that the probability that a particle does not return to the condensate 
site it has just vacated is 0(1) for all walks except for the symmetric random walk in 
two dimensions. For symmetric random walks in two dimensions there are logarithmic 
corrections: the escape probability is 0[(lnL) _1 ], where L is now the number of sites in 
the system. Therefore the upper critical dimension for the coarsening of the ZRP is two 
and above two dimensions, we expect the coarsening to be determined by the mean-field 
results for the scaling function and exponent, 5 = 1/2. 

The coarsening time scales can be summarised as follows 



for asymmetric dynamics. 

Note that mean-field analysis reproduces the correct exponent for the asymmetric 
one-dimensional case. The reason for this is that random walk arguments yield the 
same result, 5 = 1/2, whenever the escape probability is of order 0(1), and this is the 
situation in the one dimensional case with asymmetric dynamics. 

5. Heterogeneous system 

5.1. Condensation 

We now to turn to mechanisms of condensation and the associated coarsening dynamics 
in the heterogeneous ZRP. In a heterogeneous system the single-site weights fi(n) depend 
on the site I. Taking the most general case discussed in section 12*31 f where the hop rate 
from I to k is ui{ni)Wki) we have 




L 3 in 
L 2 hiL in 
L 2 in 



d= 1 , 
d = 2 , 
d > 2 , 



for symmetric dynamics, and 
r ~ L 2 Vd , 



(71) 



n 



(72) 



28 



Working in the grand canonical ensemble as in (J4.1j) one finds 



L 



Z L (z)= £ z^l[f l (m l ) = l[F l (z) (73) 

{m ; =0} 1=1 I 



where 



F l (z) = J2^ Mm), (74) 

m=0 

and (J52|l becomes 

""l^W (75) 

To simplify things we consider a one-dimensional lattice so that si = 1 and we also 
take the hop rate of the form ui(n) = Ui i.e. there is no dependence on the number of 
particles at the departure site. In this case is given by 

/,(») = (1)" . (76) 

At this point the mapping to an ideal Bose gas is evident: the N particles of the ZRP 
are viewed as Bosons which may reside in L states with energies E\ determined by the 
site hopping rates: exp(— (3Ei) = 1/ui. Thus the ground state corresponds to the site 
with the lowest hopping rate u min . 

We can sum the geometric series (|74*jl to obtain Fi = 1/(1 — z/u{) and 

>-r?^b (77 > 

Thus the maximum allowed value of z is u m { n . In the large L limit we may proceed as 
in the usual theory of Bose condensation [01] and write 

= !?z*eL + r duV{u) (78) 

L /,, . u — z 

where (n min ) is the average number of particles at the slowest site and V{u) is the 
probability distribution of site hopping rates. Interpreting V{u) as a density of states, 
equation (|78j) corresponds to the condition that in the grand canonical ensemble of an 
ideal Bose gas the number of Bosons per state is p. If the density of states vanishes as 

U ^\ ^min 

V{u) ~ (u - u min y 7 >0 (79) 

then the second term in (j7Sjl will be finite asz/ u m i n and will define p c . For p > p c , 
the first term in f)78j) must be finite thus (n m ; n ) ~ O(L). This is the mechanism of 
condensation underlying heterogeneous models of networks, for example [TTj . 
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5.2. Single defect Site 

A very simple example, which serves to illustrate the mechanism of condensation, is 
to have just one 'slow site' i.e. u% = p < 1 while the other L — 1 sites have hopping 
rates u\ = 1 when / > 1 6J. The simplicity of this system actually allows the canonical 
partition function (j^J to be calculated. This is easy to write down as 

N L 

Z l>n = J2p~ H E + (80) 

n=0 n2--.Ti.Xj fi=2 

where n counts the number of particles on the slow site. Identifying the sum over 
ri2 . . as the number of ways of distributing N — n indistinguishable balls in L — 1 
boxes implies 



n=0 v — 

In order to determine which terms dominate the sum (j81|) we look for the stationary 
point of the summand, which occurs when the ratio of consecutive terms n,n — l is equal 
to one: 

+ . = 1 . (82) 
p(N + L-n-l) K ' 

For L large, we find a solution n ~ L(p — p c ) when p > p c = p/(l — p). For p < p c , 

however, we do not find a solution for positive n therefore the maximum is the boundary 

term n = and the sum will be dominated by terms n ~ 0(1). Thus p c = p/ (1— p) is the 

critical density above which we have a condensate: we see clearly that in the condensed 

phase the defect site serves to absorb the excess number of particles L(p — p c ); in the 

fluid phase the number of particles at the defect site is 0(1). 

In the fluid phase we evaluate the asymptotic behaviour of the partition function 

by expanding for small n, since these terms dominate the sum: 

JV 



L J(N + L)^ P \ 



N+L ^ »' (83) 



L ) (l + p)(p c -p) ' 



In the condensed phase one approximates the sum by an integral and uses Stirling's 
formula to obtain 

Z^^L* f^L(I±^Y /2 (84) 
Jo V2n \ p-x J 

x exp L [(1 + p — x) ln(l + p — x) — (p — x) ln(p — x) — x lnp] 

where x = n/L. Then the integral is dominated by x ~ p — p c and expanding to second 
order and performing the Gaussian integral yields 

^^(l-p)- (i V. (85) 
From expressions (|53*l and , the mean hop rate in the two phases is given by 

(t») = rf- P<p/(1-P), (86) 
1 + p 

(it) = p for p>p/(l— p). (87) 
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These expressions for (u) are easily understood as the hop rate being limited by the 
condensate in the condensed phase and by the probability of empty sites, p(0) = 
1/(1 + p), in the fluid phase. 



5.2.1. Defect with n- dependent hop rate Recently the case of a single defect site with 
an occupation number dependent hop rate u\(n) has been considered. This allows one 
to investigate the interplay between condensation driven by particle number dependent 
hop rates and condensation driven by heterogeneity [65J. 

Again one can write down the canonical partition function. Taking 

( b\ 
uAn) = p [ 1 H — for n > 
V nj 

with p < 1 and u\ (n) = 1 for I ^ 1 yields 

f _ n\ fN + L-n- 
z ^~\ v TfTWX L-2 

71=0 

where (1 + b) n is a Pochhammer symbol defined in (|57|) . One can again determine the 
stationary point of the summand by setting the ratio of consecutive terms equal to one, 
resulting in a quadratic expression 

n (N-n + 1) 



p(b + n) (N + L-n-l) 
A real, nonnegative solution of this quadratic exists for 



9) 



P > P2 = + 7^- [% + °^I L ) ( 9 °) 



1 — p 1 — p \L 

As L — ► oo one has a condensation transition at p c = pj (1—p), which is independent of b. 
Thus the condensation is triggered by the site being slow rather than the n dependence 
of Ui(n). However on any large but finite system (with b > 0) the condensed phase 
does not appear until p 2 = p c + 0(L~ 1 / 2 ). Thus on a finite system the fluid phase 
continues above the thermodynamic critical density. This also results in the mean hop 
rate going above its predicted thermodynamic value. At p = pi one sees the system 
wander between a fluid state and a condensed state 



5.3. Dynamics 

The coarsening behaviour of the heterogeneous ZRP can be obtained by using a scaling 
argument to predict the growth law for the typical condensate size. We consider the 
case where the hop rates are drawn from the distribution with support u G [u m i n , 1] 

V ^ = n .A +7 w + > " Mmin)7 ' 7 > ° • (91) 

The argument is rather similar to the homogeneous case. After an initial nucleation 
regime, most of the system relaxes to the steady state except for a finite number of sites 
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which contain the excess particles in the system. Eventually, only two condensate sites 
will remain and these will be the sites with the two slowest hop rates. These two sites 
will be separated by a distance of order 0{L) (if they were separated by a distance of 
order 0(1) they would effectively act as a single slow site). This implies that the second 
slowest site must contain a number of particles of order O(L), since it is the slowest site 
encountered by particles in a finite fraction of the system. This is to be compared with 
the steady state, in which the j-th slowest site contains a number of particles of order 
(9((L/j) 1 ^ 1+7 ' ) ) which is subextensive for 7 > 0. This scaling can be inferred from an 
extremal statistics argument. 

In the case of asymmetric dynamics fZHl EH ESJ , the second slowest site gains 
particles from the left with rate u m { n , but loses particles to the right with rate U\ < M m i n , 
thus it loses particles with a net rate Am = u\ — u m i n ~ 0(L _1 ^ 1+7 ^). Therefore the 
time t this site takes to lose all its excess particles (i.e. excess above its typical steady 
state occupancy) scales as 0(L/Au). Thus 

T ~ L^'^ , (92) 



(93) 



and the normalised mean condensate size grows like 

1+7 

(m(t)) ( t\ 2 +-> 
(p~Pc)L \t) 

so the exponent 5 = (1 + 7)/(2 + 7) for asymmetric dynamics. 

For symmetric dynamics a similar argument can be used except that now the 
coarsening timescale r is augmented by a factor L due to the probability that a particle 
returns to the condensate site it has just vacated rather that reaching the slowest site 
[IE]. Thus 

T ~ L(3+27)/(l+7) ; (94) 



and the normalised mean condensate size grows like 

(95) 



1+7 

(m(t)} ft\ 3 + 2 ~< 



(p-p c )L 

so the exponent 5 = (l + 7)/(3 + 27) for symmetric dynamics. This was shown for 7 > 1 
in J7B], however it in fact holds for 7 > [77j . 



6. Generalisations 

6.1. Two-species ZRP 

In view of the recent interest in two component models [78 a it is natural to generalise 
the ZRP to a model with two species of conserved particles. In this way, one can explore 
the role of conservation laws in the ZRP. In particular, considering the generic nature 
of the condensation mechanism in the single species ZRP, it is of interest to look for 
other mechanisms of condensate formation which may have some generic applicability. 
Moreover, as discussed in section El the behaviour of shaken granular gases in which the 
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grains come in one of two sizes (i.e. a bidisperse system), and the behaviour of networks 
with directed edges, can be understood in terms of a two species ZRP. The particle 
dynamics we consider are chosen such that the evolution of the two particle species are 
coupled, and the nature of the coupling can be chosen such that condensation of one 
of the species is induced by the other. Such an interplay arises in models of particles 
moving on an evolving disordered background [73 EH EJ] • We discuss how this interplay 
leads to new mechanisms of condensation and rich coarsening behaviour. 

In one dimension, we define the two-species ZRP on a lattice containing L sites and 
with periodic boundary conditions. At each site /, there are n t particles of species A 
and mi particles of species B. The total number of species A particles in the system is 
N and the total number of species B particles is M. These particles hop from site I to 
the neighbouring site to the right, species A with rate u(ni, mi) and species B with rate 
v(n h mi). 

6.1.1. Steady state factorisation condition Happily, one of the most useful features 
of the single-species model — the factorised steady state — remains a feature of the 
two-species model provided the dynamics satisfy a certain constraint (HSUHSj- That is, 
one can still express the probabilities P({ni}; {mi}) in the factorised form: 

L 

P{{m}- {mi}) = Z^ NiM n f(m, mi) , (96) 

i=i 

provided the hop rates satisfy the constraint 

u(m,mi) _ v(ni,mi) 
u{n h mi-l) v(rii-l,mi) ' 

for ni,mi ^ 0. The choices of u(ni, 0) and v(0,mi) remain unconstrained. When the 
hop rates satisfy (JH7J, we can find single-site weights f(n, m) given by 

n m 

/(7i,m) = IjKi,m)]- 1 IJ[i;(0,j)]" 1 , (98) 

i=l j=l 

and with /(0, 0) = 1. We remark that the hop rates do indeed play a symmetric role in 
(|9~5j) . but that this symmetry is obscured within the constraint equation (f9*7[). Finally, 
the normalisation Zl^n,m, equivalent to the canonical partition function, is given by 

Zl,n,m = ^2 H/( ra '' m « )5 IJ^ni - N \ 5 - M ) > (99) 

{ni},{ mi } 1=1 \l=l J \l=l J 

where the 5-functions ensure that we only sum over configurations with iV particles of 
species A and M of species B. 

The steady state (jHEJ) to (jHSJ) for this model can be derived in much the same way 
as in the single species case. One begins by writing down the steady state condition on 
the probabilities P({ni}; {mi}). This has the form 

L 

= 2J [{u(7i i _i + l,m J _i)P(. • • ,71^ + 1, ni-1, . ,m;_i,m z , . . .) 
1=1 
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-u{n h mi)P(. . . , n/_i, n h , m { _i, m h . . .)} 0(nj) 

+ {u(n i _i,m i _i+l)P(. . .,ni-i,ni, . , m z _i + l, m/-l, . . .) 

-■u(nj, mj)P(. . . , rcj_i, n/, , m/_i, raj, . . .)} 0(raj)] (100) 

where the Heaviside function 0(xj) ensures that site / is occupied in order for a particle 
either to have arrived or to be able to vacate there. Now, the first term on the rhs 
of (jlOO)) is a gain term due to an A particle hopping into site I from site l — l and the 
second term is a loss term due to an A particle hopping out of site I; the third and fourth 
terms represent analogous processes for the B particles. As in the single-species case, 
we simply insert the factorised form (J35j) into (jlOOj) . but now we ask that the gain and 
loss terms due to the dynamics of the A particles cancel independently of the gain and 
loss terms due to the dynamics of the B particles. We look to achieve this cancellation 
for each term I in the sum separately hence 

u(ni, mi)/(n,_i, mi-i)f(n h mi) = u(nj_i+l, raj_i)/(nj_i+l, mj_i)/(n J -l, mi) , (101) 
for all ni ^ 0, and 

v{n h mi)f{ni^rni-x)f{n h mi) = u(ni_i,m I _i+l)/(n l _i,m { _i+l)/(n I ,m l -l) , (102) 

for all mi ^ 0, for all values of I and after cancelling common factors. These equations 
in turn imply that 

u(n h mi)f{n h mi) u(nj_i+l, mj_i)/(Tu_i+l, m z _i 



f(m-l,mi) /(m-jjmi-i) 
v(ni,mi)f(ni,mi) u(nj_i,mj_ 1 +l)/(nj_i,mi_ 1 +l) 



constant , (103) 
constant , (104) 



f(n h mi-l) f(m- 1 ,m l ^ 1 ) 
for all I. Both constants are set equal to unity without loss of generality, since they only 
appear as overall factors in the normalisation. The two relations (|103|) and (|104|) are 
iterated to yield (|§H|). However, (|103j) and ()104|) also imply the constraint (|57j) on the 
choice of u(ni,mi) and v(ni,mi). To obtain this constraint, we use ()103j) and ()104j) to 
obtain two expressions for f(ni,mi) in terms of f(ni — 1, mi — 1): 

,/ \ /(ni-l,m/-l) /(n,-l,mj-l) , . 

f{n h mi) = — — r - = — — r. (105) 

u{n h mi)v{ni-l,mi) u{n h mi-l)v[n h mi) 

Since both of these expressions must give the same result the hop rates are required 
to obey the constraint (|97|). Equation ()105J) also indicates how we might interpret 
physically the constraint (|37j) : it suggests that the relationship between /(nj — 1, mj — 1) 
and /(nj, m/) is the same when an A particle hops from I and then a P, as when the B 
particle hops before the A. In other words the single-site weight f(ni, mi) is independent 
of the order in which the particle species arrived at I. 

As in the single-species case, it is straightforward to generalise the derivation of 
the steady state to an arbitrary lattice, including generalisations to disorder in the hop 
rates, partial asymmetry in the dynamics and any dimension. It is also possible to 
generalise to any number of species, Q say, in which case each species satisfies Q — 1 
constraints of the form (|57j) , one for every pair of species 
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6.1.2. Condensation We now consider how the interaction of the two species may lead 
to a condensation transition as a function of the particle densities, 

N M , , 

Pa = -j- and p B = — , (106) 

for species A and B respectively. We consider a case where the dynamics of one of the 
particle species — the B particles — depends only on the number of particles of the 
other species at the departure site. Hence we take 

v(n,m) = 1+ ^ , (107) 
[n + l) 7 

for all n and for m > 0, where c, 7 > 0. This choice for v(n, m) now determines the 
m-dependence we must take for u(n,m) through the constraint (|97j): 

/l + c \ m 

u(n,m)={ l + T ) ■ ( 108 ) 

We could choose to multiply u(n,m) by some function of n and still satisfy (}9"T|) but 
we will proceed with just the rate f)108j) for simplicity. The single-site weights f(n,m), 
determined by (jHEJ), are given by 



As in the single-species case, it is simplest to work within the grand canonical 
ensemble in order to demonstrate condensation. The grand canonical partition function 
is defined as 



00 00 



Z L ( Z , V)=J2Y1 Zn y mZ L,n,m , (HO) 

71=0 771=0 

where we now have two fugacities, z and y, which fix the particle densities through 
zd\nZ L (z,y) y d\nZ L (z,y) 

PA = L—^z ' PB = L dy ■ (1U) 

Now, we use ([99)1 and perform the sums over n and m in (jllOj) to obtain 

L 

Z L (z,y)= ^ ini y Limi X{f{n u m l ) = [F{z,y)} L , (112) 

{ni},{mi} 1=1 

where 

00 00 

F^l/j^^.y/^m), (113) 

71=0 771=0 

and the densities are determined by 

d\nF(z,y) d\nF{z,y) 

pa = z—^- , pb = y^^~ ■ ( 114 ) 

In the single species case, the critical density was determined by the properties of F(z) 
and its derivative when z assumed its maximum value (3. A similar picture emerges 
here, which we illustrate for / '(n,m) given by (|109|) . 
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Figure 6. The solid lines represent schematic dependences of the particle densities 
Pa and ps for fixed values of z and as a function of y. The dashed line in the right 
hand graph illustrates how ps varies as a function of z and y given that pa is fixed. 



Using (Jl(J9j) in (J113|) . and performing the sum over m, one obtains 
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F(*,y) = Y>" {1 tZ + " , (H5) 
^ (l + np(l-y)+c } 1 ; 

d , x A „ (l + n) 7 + c . . 

^ ^ (l + n)7(l- y ) + c ' 1 > 

vlnz, y) = E(i + n)^ ,/ ((1 + " )7 + c) 12 ■ (H7) 

The radii of convergence of these series are z — 1 and y = 1. We can now use equations 
fl 11 5 j) to (|117j) in (|114|) to find z and £/, given values for pa and p#. In figure El we 
illustrate how pa and ps depend on y for fixed z. This figure is obtained by noting that 
for a given z, both pa and p# are monotonically increasing functions of y. Then we can 
infer the behaviour shown in figure El by analysing (|114|) in the four limits y — > 0, z — > 0, 
y — ► 1, and z — > 1. The existence of a condensation transition is demonstrated by 
considering a system containing a density pA of A particles, represented by the dashed 
lines in figure |H1 The solution of ()114|) requires that z and y lie in the range Z\ < z < z 2 
and < y < 1. In this range, p# increases monotonically from p# = 0, where z = z 2 
and y — 0, to its maximum value, where y = 1 and 2 = Z\. Thus, for every finite value 
of pa there exists a finite maximum value of ps above which we can no longer solve 
flll4jh When ps exceeds this maximum (|114jl cannot be satisfied and a condensation 
transition ensues. Hence we obtain the phase diagram shown in figure For 7 = 1 the 
critical line is given simply by 

Pb = (1 + Pa)/c (118) 

Thus the model exhibits a phase transition between a fluid phase, in which both 
particle species are distributed exponentially, and a condensed phase, in which p B 
exceeds a maximum critical value and the excess B particles condense onto a single 
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Figure 7. Phase diagram for the two-species zero-range process with 7 = 1. 

site. The number of B particles in this condensate is proportional to the system size 
L — in the limit N, M, L —>■ oo it contains a finite fraction of an infinite number of 
particles. However, we note that for m — > oo, u(n,m) — > if n is finite. But this is 
not consistent with continuity since the current of A particles away from the condensate 
site is finite therefore u(n, m) at the condensate site must also be finite. Thus it must 
be that n — > oo at the condensate in order that the current out of the condensate site 
is finite. For large n and m, u(n,m) ~ exp(— r ym/n 1+1 ) and since m is a number of 
order m ~ O(L) at the condensate site, n must be of order n ~ O^L 1 ^ 1 ^) for finite 
u(n,m). Thus the condensate of B particles is sustained by a subextensive condensate 
of A particles. 

The conditions on u(n,m) and v(n,m), and therefore f(n,m) also, leading to 
condensation have been considered for general rates in |HJJ. Additional condensation 
transitions can be obtained by supplementing the A particle hop rate ()108|) with a factor 
(1 + b/n). In this case two additional condensed phases arise, one in which only the A 
particles condense, and one in which both particle species condense — thus this phase 
also exhibits induced condensation of one species by the other. 

6.1.3. Dynamics of condensation The coarsening dynamics in the condensed phase 
shown in figure can be obtained by adapting the random walk arguments given for 
the single species case [So] . 

Thus, as before, we consider a system most of which has relaxed to the steady state, 
except for a finite number of condensate sites. The condensate sites contain a number 
of A particles of order n ~ (9(L 1 ^ 1+7 - ) ) and a number of B particles of order m ~ O(L), 
i.e. m ~ n 1+1 . Therefore, from (jl07|) . the B particle hop rate decays to its steady 
state value as c/m 7// ^ 1+7 \ so the typical timescale over which B particles escape from 
condensate sites scales as L 7 ^ 1+7 \ Then the time it takes a condensate site to lose all of 
its B particles scales as L 1+7//< - 1+7 \ This sets the coarsening timescale, thus the growth 



37 



of the mean B particle condensate size follows the scaling law (|67|) with exponent 

j =r^- < 119 > 

A prediction for the exponent 5 can also be obtained in the case that the dynamics 
are symmetric. The dynamics are altered in the same way as in the single species 
case, in that only every L-th particle on average reaches the next condensate site rather 
than returning to the one it has just departed, thus the typical timescale over which 
B particles escape from condensate sites is a factor of L longer: it scales as L 1+7 /( 1+7 ). 
Therefore the time it takes a condensate site to lose all of its B particles scales as 
^2+7/(1+7) anc j coarsening exponent is given by 

1+7 

5 = — V- . (120) 
2 + 3 7 v ; 

It is interesting to note that these exponents are the same as those obtained for the 
heterogeneous ZRP with hop rates drawn from the distribution (c.f. section EHH) 

v{u) = [(7- 1 + i)/(i - u min y- 1+1 \(u - u^y 1 . (121) 

To consider why this might be the case, one views the dynamics defined in fll07H08|) 
as a model of particles (the B particles) moving on an evolving disordered background 
(given by the A particles). By the time the coarsening regime has been reached, at 
the condensate sites the evolving disorder is effectively quenched. Therefore it is not 
necessarily surprising that the two models exhibit similar coarsening behaviour for some 
distribution V(u). The reason the form ()121)) is the relevant one for the rates (jl07|108j) 
is as follows. In the heterogeneous system, as discussed in section 15. 3| the coarsening is 
governed by the difference between the hop rates at the two slowest sites in the system, 
An. For the distribution (j!21j) . Am ~ £^ 1 /( 1 +t 1 ), This rate separation plays the role 
of the background of A particles in the two-species model. The remaining contributions 
to the coarsening time scale in the heterogeneous model are then determined by the 
symmetry of the hopping dynamics, i.e. the coarsening time scale is given by a factor 
of order L for asymmetric dynamics, or a factor of order L 2 for symmetric dynamics, 
multiplied by the inverse rate separation. This leads to the same exponents as those 
obtained in the condensed phase of the two-species model. 



6.2. Urn models and the Misanthrope process 

A natural extension of the ZRP dynamics is to hop rates which depend not only on 
the occupancy of the departure site, but also on that of the target site. For example, a 
recent application of such models is to rewiring dynamics in networks as we discuss in 
section 16.2.41 Models of this kind are often referred to as 'urn models'. As reviewed in 
[HE] urn models comprise balls distributed amongst a number of boxes with conserving 
dynamics for the exchange of balls. These dynamics are usually defined on a fully- 
connected geometry and obey detailed balance with respect to some energy function 
defined as a sum of single-site energies. 
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On the other hand the Misanthrope process, introduced in the mathematical 
literature in j^Zj, is a model with dynamics defined without reference to an energy 
function. In one dimension its dynamics are defined in exactly the same way as the 
ZRP except that now the hop rate depends on the occupancy of both departure and 
target sites. Thus a particle hops from site I to site Z + l with a rate u(ni, n/+i). 

6.2.1. Steady state factorisation condition It turns out that the steady state still has 
a factorised form, provided the hop rates satisfy a constraint. To see this, consider the 
condition on the steady state probabilities P({ni}): 

L 

= J2 [tt(n,_ 1 +l,n,-l)P(---,n,_ 1 +l,Tu-l,--O0(n I ) 
1=1 

-uim^n^PHm^eim-t)) . (122) 

To look for a factorised steady state of the form (J2J) we introduce a counterterm which 
cancels under the sum and then look to cancel each term in the sum separately, hence 

/(n,_i)/(n,) - fin^fim) = ufa-i+l, n,-l)/ (n w +l)/(n,-l)0(rO 

-u(ni-i,ni)f(ni-i)f(ni)9(ni-i) , (123) 

having cancelled common factors (a product over the functions f(n,k) at all sites k ^ l—l, 
I), where f(n) is some function to be determined. The lhs of (|123|) represents the 
counterterm. We remark that this cancellation mechanism is identical to that of the 
so-called 'Matrix Product Ansatz' in the case where matrices are replaced by the 
functions /(n) and f(n) |HBj- 

In order to satisfy ()123|) for all values of n\_\ and n\ there are now three cases to 
consider: 

In this case 

/(0)/(n,)-/(0)/(n,)=tt(l,n,-l)/(l)/(n,-l) , (124) 

• ni-i 0, n t = 
In this case 

/(ni_i)/(0) - /(n,_i)/(0) = -u{rn- u 0)/(n,_i)/(0) . (125) 

The lhs of fll24j) and (|125|) are both of the same form, therefore we are able to 
eliminate f(n) from these two equations which yields 

u(n,0) /(0) 

for n > 0. This recursion is easily iterated to obtain the factors f(n) in terms of 
the hop rates: 
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a ) _u(3,l) 

4i(l,0) 



o 

_ o_ o_ 

Site: 1 2 3 4 5 



u(3,l) u(l,0) 



b) 



Particle: 1 2 3 4 5 

Figure 8. Mapping between the Misanthrope process and the asymmetric exclusion 
process 



We can now substitute ()126j) into ()124|) or (jl25j) to find the expression for f(n), 
hence 

/(n) = /(n)[/(0)-tt(n,0)]. (128) 
• ni-i ^0,ni^0 

We have used the first two cases to derive expressions for f(n) and f(n). Thus we 
must now satisfy flUS)) under substitution of QMfy and {HHP - We find that flUHD 
is indeed satisfied provided the hop rates satisfy the condition 

u(n, m) - u (n+l,m- 1) , <^ n H m ^) = u , 0) _ u , 0) _ ( 129 ) 

u(n + 1, 0)u(l, m — 1) 

This condition reduces to two separate conditions 

u(l,n)u(m,0) 

u{m,n) =u(n + l,m-l)— V - ' (130) 

u(n + 1, 0)w(l, m — 1) 

u(n, m) — u(m, n) = u(n, 0) — u(m, 0) (131) 

for all values of n and for m > 0. 

Hence the steady state of the Misanthrope process is given by the factorised form (0) 
with f(n) given by (|127|) provided the hop rates satisfy the constraint ()129|) . We remark 
that if we take u(n, m) to be a function of n, the number of particles at the departure site 
only, we recover the ZRP, albeit via a slightly different derivation to the one presented 
in section El 



6.2.2. Mapping to the KLS model Just as for the ZRP, the Misanthrope process can be 
mapped onto an exclusion process. This mapping is illustrated in figure El It is achieved 
in the same way as for the ZRP — sites in the Misanthrope process become occupied 
sites in the exclusion process and particles in the Misanthrope process become vacant 
sites in the exclusion process — the only difference is that now the hop rates depend 
not only on the distance to the next particle in front, but also on the distance to the 
previous particle behind. 



40 



An example of this mapping is between the Misanthrope process and the KLS model 
The KLS model is an exclusion process in which particles hop from site i to site 
i — 1 with a rate depending on the occupancy of sites i — 2 and i + 1. The dynamics are 
defined by the processes 



1011 1101, 1010 1100 

a 3 



0010 0100, 0011 -^0101, (132) 

where we use a '1' to represent a particle and a '0' to represent a vacancy. The steady 
state probabilities P({r)i}), where the occupation variables rji = or 1, can be written 
in the form 

1 / L+N \ 

(KLS) eX P 
^L,N \ i=l / 

for a system with L particles and N vacancies, where Z^^f * s a normalisation. This 
steady state holds provided [HU] 

«2 = e^Q?4 , «2 + «4 = Oil + «3 . (134) 

The mapping to the Misanthrope process is achieved by making the identification: 
u(n,m) — a 3 for n > 1, m > 0; u(n, 0) = ct 4 for n > 1; u(l, m) = a 2 for m > 0; 
m(1,0) = «i. With this choice of u(n,m), condition ()129j) reduces to ()134|) . This leads 
to a Misanthrope process for N particles on L sites with the steady state (J2J where /(n) 
is given by 

/(*) = (y|y) /(°) e/3(n_1) for » > 1 » ( 135 ) 



and where substitution of the Misanthrope rates into ()129|) yields the second of the 
conditions in (J 134)) . Thus, under the mapping the Ising-like form of the steady state 
()133|) changes into a simple factorised form. 

6.2.3. Misanthrope process on fully- connected geometry In this subsection we consider 
a model where a particle can hop to any other site with a rate which depends on 
the occupations of the departure and destination sites. Thus we are considering a 
Misanthrope process on a fully-connected geometry. 

Since the sites are fully connected the steady state weight should be invariant under 
all permutations of the sites. Note, however that on a finite system this does not imply 
a factorised steady state. This is perhaps counterintuitive as one might expect the fully- 
connected geometry to yield 'mean-field' results and the steady state to factorise. On 
the other hand, in the limit of infinite system size mean field results should be recovered 
and it would be interesting to demonstrate that this is the case. 

On a finite system, it turns out that the steady state has a factorised form only if 
the hop rates satisfy a constraint. This constraint is less restrictive than the asymmetric 
one-dimensional case discussed above. When this constraint is satisfied detailed balance 
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is obeyed. Thus, in this case the Misanthrope process is equivalent to the urn models 
reviewed in [86 . 

To see this, consider the condition on the steady state probabilities: 

= K n ^+ 1 ' n fc _1 )- P (' ' ■i n i + 1 ' • ■ ■ ,n fc -l, ■ • •) - u(n h n k )P({ni}) 

+ u(n k + l,ni-l)P(- ■ - jTij-l, • ■ - ,n fe + l, ■ ■ •) -u(n k , n z )P({n z })(|136) 

where the sum runs over all pairs of sites I, k. The assumption of a factorised steady 
state yields 



l^k 



' /( Wi + l)/ (7lfc -l) 

r*(n 1+ l,n»-l) /(n;)/K) -^n,) 



/(n fc + l)/(n ; -l) 

+w(n fc +l,n z -l) — — u(n k ,ni) 

f{n k )f(ni) 



(137) 



It turns out that to satisfy this one must have detailed balance 

u(ni+l,n k -l)f(ni+l)f(n k -l) = u(n k ,ni)f(ni)f(n k ) (138) 

which yields 

/(»)=/(») f^y f[ ^^^r < 139) 

\/(°)/ t=l u ( m '°) 

provided the following condition holds for n > 1, m > 

w(n, 0)u(l,m) 

m) = U (l,„-lMm + l,0) u(m + hn ~ 1} • (M0) 
which is one of the two conditions (jl3())l required for the one-dimensional asymmetric 
case to factorise. A simple way to guarantee the equality of ()140|) is for u{n,m) to 
factorise into a function of n and a function of m. This is the case of interest for 
rewiring networks discussed in the next subsection. 



6.2.4- Mapping to network dynamics Recall that a network is a collection of nodes 
connected by edges with dynamics for rewiring of edges (see section 13. 2j) . Here we 
discuss how these dynamics correspond to Misanthrope dynamics. 

Let us fix the number of nodes in the network to be L and the number of 
(undirected) edges to be N. A general rewiring dynamics is as follows 39]: an edge is 
chosen at random and one end of the edge initially connected to node k is rewired to 
another node k! with rate (p(k')h(k). An increasing function 4>(k') describes preferential 
attachment whereas a decreasing function h(k) describes enhanced detachment of an 
edge from nodes with low degree. (In the network literature f(k) is usually used to 
denote <fi{k) but we have chosen to use <f>(k) to avoid a clash of notation with the single 
site weights f(n).) Thus the overall rate of a rewiring event such that two nodes with 
degrees k, k' end up with degrees k — 1, k' + 1 is proportional to k(f>(k')h(k). Note that 
the factor k stems from choosing an edge at random. The dynamics clearly preserves 
the number of edges and thus generates the 'canonical ensemble'. Other ensembles of 
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graphs are discussed in j!2] and jlSj- Note, however, that the dynamics may generate 
self-connected nodes (tadpoles) and multiple edges between two nodes (melons). It has 
been argued that these graphs have little effect on the resulting ensemble |3Hj . 

We now think of the number of edges connected to a node as the number of particles 
at a site in the Misanthrope process. The rate at which a particle hops from a site with 
k' to a site with k particles is then u(k, k') = kh(k)4>(k') . Also note that an edge may 
be rewired from one node to any other, thus we are considering Misanthrope dynamics 
on a fully-connected geometry. We may therefore invoke the results of 1(1 2. HI which tell 
us that if u(k, k') factorises, which is the case here, we have detailed balance and the 
steady state factorises. The single-site weight becomes 



Recalling the results of section l4~2l we see that for condensation to occur we require 
that mh{m) / 4>{m — 1) decay more slowly than 1 + 2/m. Setting h(k) = 1 gives the 
preferential attachment case. Thus the preferential attachment function <j)(m)/m should 
increase more quickly than 1 - 1/m to have condensation. 



An alternative mapping from the dynamics of a network to a ZRP is proposed in 

EH- 



6.2.5. Relation to Backgammon model The Backgammon model, introduced in [92J, 
is an urn model which can be related to the Misanthrope process. There has been 
considerable interest in the Backgammon model as a simple model of glassy dynamics. 
Defined on a fully-connected geometry, the idea is that particles hop between sites and 
there is an energy cost associated with occupied sites. The hopping rate from site 
I to k is therefore a product of factors for each site and the existence of an energy 
function guarantees that condition (|130j) is satisfied. Since there is an energy function, 
the dynamics can be parameterised by a temperature and at very low temperatures the 
steady state will be dominated by configurations with very few occupied sites. Glassy 
dynamics are exhibited in the relaxation towards the steady state, where the elimination 
of occupied boxes becomes slower and slower. 

A simple example of Misanthrope dynamics corresponding to the Backgammon 
model is u(n,m) = 1 + (e~P — l)5 m ,o for n > 0, where j3 is the inverse temperature. 

6.3. General mass transfer model 

Given that the steady state of the ZRP factorises it is natural to ask under what 
conditions may a factorised steady state be admitted. Recent work [§311111] has answered 
this question and revealed an appealingly simple condition for factorisation within a 
framework that encompasses a wide range of models. 
The work [531 H31 generalises the ZRP in three ways 

• The mass is a continuous variable 




(141) 



43 



• Arbitrary amounts of mass may move from one site to another 

• The dynamics comprises discrete (or parallel) timesteps 

The class of models considered are defined on a one-dimensional lattice of L sites 
with periodic boundary conditions (site L + 1= site 1); associated with each site is a 
continuous mass variable m,, i — 1 . . . L. The total mass is given by M = Yli=i m i- 

The dynamics is defined as follows: at each time step, at each site i, mass \ii drawn 
from a distribution </>(/Xj|mj) 'chips off' the mass rrii, and moves to site i + 1. Thus all 
sites are updated in parallel and mass is transferred simultaneously between sites. 



6.3.1. Steady state factorisation condition The central result of [93 is a necessary 
and sufficient condition for the steady state to factorise which states that the chipping 
function must be of the form 

v (n) w (m — /i) 



[v * w\ [m) 

In this case the single-site weights become 

f(m) = [v * w] (m) , 
where in expressions f!142|) and (j!43|) the convolution is defined as 



[v * w] (m) 



dfi v (/i) w {m — n) 



(142) 



(143) 



(144) 



Equation ()142|) expresses that the hopping function should be a product of a function 
of n, the mass which moves, and a function of m — fi, the mass which remains, divided 
by a normalisation which is a function of m. 

The proof that ()142j) is necessary and sufficient is presented in Here we content 
ourselves with checking that it is a sufficient condition. 

The steady state probabilities satisfy the following condition 

poo /-mj L 

P(m) = JJ / dm[ I dfii <p(fii\m^Y[S{m k - m' k + fi k - fi^) P(rr^) . (145) 
i=i Jo Jo k=i 

The r.h.s. of this equation is to be understood as integrating over all possible 
configurations before the update: site / contained mass m[ = mi + fii — and mass 
111 moved to site I + 1 at the update. Assuming a factorised steady state, equation (|145|) 
reduces to 

L L 



11/1'"/) = II 



i=i 



[ J 4>{^k\m k + fi k - ii k -i)f{m k + n k - /i fe _i) 



Inserting (|142|143|) the r.h.s becomes 



n 

i=i 



d/i; 



Y[v(fx k )w(m k - = Y[ [v * w) (mi) 



(146) 



(147) 



i=i 



where we relabelled k — I + 1. Thus ()145|) is satisfied when ()142I143|) hold. 
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6.3.2. Continuous time limit The continuous time case is a limit of the discrete time 
case when we let the time step be dt and set v (//) = <5(/i) + x(/^)dt then take the limit 
dt — > 0. Condition (jl42j) reduces to requiring hopping rates 7(/x|m) of the form 

. . x{n)w(m — y) ,,. a . 

7(^M = 7 \ 148 

u;(to) 

in which case /(to) = w(m). 

It is useful to check explicitly that this works in a way similar to section 12.21 . The 
steady state probabilities satisfy the following condition which generalises (J3J) 

L r rmi 

= ^ / d/i/_i7(/i/_!|^_i + m,_ x )P({. . . ,1711-! + Vi-urn - /J,_i, . . .}) 



Z=l 



mi 

d^7(/i;|m i )P({nJ) 



o 



(149) 



Assuming a factorised steady state and equating terms I in the sum, yields 

mi 

d//j_i7(juj_i|/^_i + mi-i)f(mi-i + /ii-i)f(mi - = 
dna(lj,i\mi)f(mi-i)f(mi) . (150) 



o 



o 

Inserting ()148j) and /(to) = w(m) yields the required identity. 

The case of discrete masses is easily obtained by restricting /i, m — [A to be positive 
integers in (j!42)l and defining the convolution as 

m 

[v * w] (m) = v (//) w (m — //) . (151) 

jU=0 

Heterogeneous case We also note that the condition (|144|) can be generalised to 
the heterogeneous case where 0,(/x|m) depends on the site /. A necessary and sufficient 
condition for a factorised steady state is that 
, , . . v (u)wi (m — u) 

0, A*M = 7 *\, , W , (152) 

where v and if; are arbitrary functions but v must be the same for each site. The 
single-site weights are given by 

fi(m) = [v * wi] (to) . (153) 

The continuous time limit implies that the site-dependent hopping rates should be 
of the form 

. . . x{a)wi{m — a) 
7 (fi to = w l j w , 154 
wi(m) 

in which case fi(m) = wi(m). 
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6.3.4- Test for factorisation Although (|142j) is appealingly simple, it does not tell us 
directly if a given hopping function yields a factorised steady state. Happily, a direct 
test is easily constructed [M]. We require that 



d_ 



d_ 



ln0(/x|/x 



cr 



h{jJ, + a) 



(155) 



i.e. the result of differentiating In <fi with respect to fi, the mass that moves, and a, the 
mass that remains, is a function of m = a + fi alone. Then from ()142|) we would have 
h(m) = —d 2 ha.f(m)/dm 2 and integrating twice yields 



/ (m) = exp 



dm' 



dm"h (m") 



(156) 



In the same way we can construct a test for the integer mass case: we require 
</>(// + l\n + 2)<f>(fj,\n) 



R{n) , 



(157) 



<p(p+ l|n+ \)<f>{n\n- 1) 
i.e. the cross-ratios on the l.h.s (defined when all of the 0's are positive) are functions 
of n alone. Then from ()142j) we would have 

R(„) - f{H + 1)2 
K[n) f(n) f(n + 2) ' 

which implies the recursion 

/(n + 2) 1 f(n + l) 



158 



/(n+1) R(n) f(n) 
Iterating ()159j) twice yields 

/fmN nn-2 

/(") 



/(o)^Vn[fi— 



for n > 2 . 



(159) 



(160) 



6.3.5. ZRP with parallel dynamics As an illustration of the use of these tests and the 
construction of factorised steady let us first consider the ZRP with parallel dynamics. 
At each time step a particle hops forward from a site with n particles with probability 
u(n). Then <p(0\n) = 1 — u(n), <p(l\n) = u(n), <ft(k\n) = for k > 1. 
The cross ratio in ()157|) is only defined for fi = 0, therefore 
u(n + 2)u(n) 



R{n) 



(161) 



u{n + l)u(n — 1) 

is automatically a function of n alone and factorisation is guaranteed. Expression (jl60J) 
becomes 



fin) = /(0) 



/(o) 



n 

n 1 



«0') 



(162) 



This expression was first derived in the exclusion process context [^3] by a more 
complicated approach. Also in that work expressions for f(m) were presented for the 
case of ordered sequential dynamics where one updates either in the forward sequence 
I — 1 . . . L or the backwards sequence I = L . . .1. 
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6. 3. 6. ZRP with transfer of more than one mass unit As a more complicated example 
we consider a generalised ZRP where mass 'chunks' of size one or two can chip off at 
each time step with probabilities Ui(n) and w 2 (n) respectively. In this case we have two 
cross ratios in (| 15 7(1 defined for /i = 0, 1. Demanding that these be equal implies 

R(n) = "i(" + 2)[l- M i(n)- M2 (n)] 
v ; ui(n + l)[l-«i(n + l)-u 2 (n + l)] v ; 

u 2 (n + 2)Mn) for n > 1( (164) 



u 2 (n + 1) Mi (n + 1) 
which reduces to 

u 2 {n + 1)(1 - ui(n) - u 2 {n)) 



A for n > 1 , (165) 



Wi(n + l)iii(n) 

where A is a positive constant independent of n. This recursion fixes it 2 (n) in terms of 
A, w 2 (l) and «i(n). Then we can use (I16U|) and (|163|) to deduce the single-site weights 



f(n) = /(0) 



(my m^t a i-^oo-^q-) , * 

V/(0)j 1 - u x {n) - u 2 (n) 11 UiO") ' 1 J 



7. Condensation We now analyse the condensation phenomena arising from /(n) 
of the form ()143jl . For brevity we restrict our attention to the homogeneous case. The 
grand canonical partition function becomes 

L 



Z(z) 



dm z m [v * w] (m) 



[V(z) W(z)} L , (167) 



where 



/"OO /"OO 

V{z)= / dfiz^vfji) ; W(z)= / daz a w(a). (168) 
Jo Jo 

As in (J32I35J) the mass density p determines the fugacity z through 

P = TW + ^W (169) 

As before the condensation mechanism is as follows: since (j!69)l is an increasing function 
of z one can always find a finite solution of of equation ()169J) provided the r.h.s. diverges 
at the radius of convergence of V or W. On the other hand, a condensation transition 
can occur if zW'(z)/W(z) and zV'(z)/V(z) have finite limits as z tends to the lower of 
the radii of convergence of V(z),W(z). 

Note the condensation mechanism may be driven by either or both of v(fi), w(o~) 
according to which of V(z), W(z) has a finite radius of convergence. The condensation 
in the ZRP corresponds to condensation driven by w(a) since in the ZRP the amount 
of mass that can chip off is restricted to one unit. Condensation driven by v(fi) would 
require a slowly decaying v(/a) which would correspond to chipping off large fractions 
of mass. One would expect that under these conditions the condensate is mobile since 
events where a large fraction of mass is chipped off will move the condensate around the 
lattice. 
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Thus in addition to the two types of condensation seen in the ZRP i.e. condensation 
of particles onto the site with the slowest hopping rate in a heterogeneous system and 
condensation due to slowly decaying hop rates on a homogeneous system, we identify 
a third type of condensation due to slowly decaying size distribution of the chunks of 
mass transferred. 

6-4- ZRP with non- conservation 

So far we have considered the ZRP and generalisations where the particle number, or 
mass, and site number is conserved. In this section we explore some cases where these 
quantities are allowed to fluctuate yet the factorisation properties are preserved. These 
dynamics allows one to generate different ensembles from the canonical ensemble. As 
we shall see, this admits a variety of new transitions. 

We restrict ourselves to periodic boundary conditions although we note that an 
open boundary ZRP has recently been considered in [96J. 

6.4-1- Fluctuating particle number and generation of grand canonical ensemble Let us 
consider dynamics where in addition to the hopping rates u(n), a particle is created 
at site I with rate c{n{) and a particle is annihilated with rate a(ni). In order for the 
steady state to remain factorised and for the particle number to fluctuate, the dynamics 
for the creation and annihilation of particles must obey detailed balance with respect 
to the steady state 

p (W}) = mqm {m 

where /(n) = f{n — l)/u{n) satisfies the usual factorisation condition (jHJ). Detailed 
balance with respect to particle non-conservation at all sites implies a(n)f(n)z n = 
c(n — l)f(n — l)z n ~ l so that 

a(n)z = u{n)c{n — 1) . (171) 

A convenient solution of this equation is a(n) = u(n) and c(n) = z. Thus, with this 
choice of rates one generates in the steady state the grand canonical ensemble of section 
IP 

If c(n)/a(n) > 1 as n —>■ oo, the particle number will diverge and a steady state 
will no longer be attained. In section 14.11 a cut-off on particle number, which implies 
f(m) = for m > N, was introduced to circumvent this problem. In order to generate 
the cut-off dynamically one sets f(N+ 1) = f(N)/n where k — > oo. Then one can take 
u(N + 1) = a(N + 1) = k -> oo. 

Finally, in section l4~Tl z was taken as dependent on the total particle number z(N). 
This can be generated dynamically by letting the creation rate be c = z(N). Note that 
such a dependence of local particle creation rate on the global number of particles in 
the system implies non-local dynamics. 
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Conserving dynamics 
ZRP: 

'o 



Exclusion process: 



u(n) 



n-1 



u(n) 



n n-1 
Non-conserving dynamics 
ZRP: 



o 




wu(n) 




Exclusion process: 



wu(n) 



n n-1 

Figure 9. Illustration of the dynamics which allows the site number and the particle 
number to fluctuate with L+N fixed. The conserving dynamics, as shown at the top 
of the figure, are the usual ZRP dynamics where particles hop to the neighbouring site 
with rate u(n). The corresponding exclusion process dynamics are illustrated beneath. 
The lower part of the figure illustrates the non-conserving dynamics: simultaneously 
a particle is annihilated and an adjacent vacant site created with rate wu(n) and the 
reverse process occurs with rate one. The corresponding exclusion process dynamics 
are shown at the bottom of the figure. 



6.4-2. Fluctuating site and particle number We now consider a microscopic dynamics 
which allows the number of sites L and the number of particles N to fluctuate under 
the constraint that the total number of sites plus particles remains fixed: L + N = K. 
The dynamics is illustrated in figure |H1 along with the mapping to an exclusion process. 

It is easy to check that this choice of creation and annihilation dynamics satisfies 
detailed balance with respect to the steady state 

P({ ni . . . n L }) oc J] [wz^fim)] (172) 
i 

where, as usual, f(n) = f(n — l)/u(n). That is, the creation of new sites balances 
annihilation of sites through 

wu{n) [wz n+1 f{n)] = [wz n f{n - 1)] [wzf{0)] . (173) 

An interesting phase transition occurs under these dynamics which is most easily 
studied in a grand canonical ensemble. We define the grand canonical partition function 
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as 



where F(z) is defined as in (|34jl . The value of 2 is fixed through the condition 

K = z-^— = (L) [(n) + f] (176) 

where 

. . dh\Z 1 , 
L = = — - 177 

ow 1 — wzF(z) 

<»> - ■ (178) 

(L) is the average number of sites and (n) is the average number of particles per site. 
When K is large, ()176j) is satisfied when either (L) or (n) is large. 

Let us take as example the solvable case of section E.4.11 u(n) = l + b/n and for 
simplicity we consider only b noninteger. 

If b < 1, F(l) diverges and we can clearly find z < 1 such that 1 — wzF(z) = 
0(1/ K). Thus for b < 1 we have from (117711781) that (L) = O(K) and (n) = 0(1). 
This implies a number of sites of order K each with a finite number of particles. 

For b > 1 we have the explicit expression F(l) = b/{b — 1). Thus if w > w c = 
(b — l)/b we can again choose z to give (L) = O(K) and (n) = 0(1). For it; < w c , on 
the other hand, we have different behaviour according to whether b < 2 or b > 2 

For b < 2, -F'(l) diverges. Therefore for w < w c we can choose z 1 such that 
(n) = 0(K). In this case (L) = O(l), thus we have a finite number of sites each with 
a number of particles of order K. 

For b > 2, -F'(l) = bj(b — 1)(6 — 2), and for 10 < w c we cannot satisfy ()176|) as 
2 /* 1. Moreover as z /* 1, both (L) and (n) are finite. Therefore we must have a 
condensation transition wherein a condensate site emerges which contains fraction one 
of the particles in the large K limit. 

The mechanism for these transitions was first studied in a simple lattice model for 
DNA denaturation [HI] . Within the mapping to an exclusion process the constraint that 
K is constant corresponds to fixed number of lattice sites with particle non-conservation. 
Thus the transition is from a state with a finite density of particle to a zero density state 
which comprises a finite number of particles with large gaps between for 1 < b < 2, or 
the condensed state with one large gap and a small jam of closely spaced particles for 
b > 2. 

It has been shown how the transition may occur in an exclusion process with local 
dynamics by using two species of particles to generate effective hop rates such that 
b = 3/2 and f(n) ~ 1/n 3 / 2 , in the spirit of section EDO [37|. 
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6.4-3. Fluctuating site number only One can also consider the case where only the site 
number L is allowed to fluctuate [98] . For dynamics in which, in addition to the usual 
ZRP hop rates, vacant sites can be created with rate w and annihilated with rate 1, the 
steady state 

L 

P L ({ni...n L }) <xw L Y[f( ni ) , (179) 

i=i 

satisfies detailed balance with respect to the creation and annihilation dynamics. The L 
superscript here is intended to emphasise that this is the probability for a configuration 
of iV particles on exactly L sites. 

This model undergoes a transition which can be analysed in a way following just 
the same steps as the previous case, so we neglect the details and quote only the results. 
We consider hop rates u(n) = 1 + b/n, and we consider the three quantities N, (L) and 

<»>• 

There exists some critical value of w, w c , such that if w > w c then we can work in the 
grand canonical ensemble and find a fugacity such that (L) = O(N) and (n) = 0(1). 
However for w < w c two cases ensue. If b < 2, one can find a fugacity such that 
(n) = 0{N) and (L) = 0(1). On the other hand, if b > 2 one finds that both 
(L) = 0(1) and (n) = 0(1). Therefore in the limit iV — > oo, a single site contains a 
fraction equal to one of the N particles in the system. 

7. Summary and open questions 

In this work we have reviewed the properties of the zero-range process and a variety of 
generalisations and related models. We have focussed on the property of a factorised 
steady state which is exhibited in many of these models: the ZRP always factorises, 
whereas generalisations such as the two-species model (section Ifi.l]) and Misanthrope 
process fsection ltT^ . where the hop rate also depends on the target site, only factorise 
when the hop rates satisfy certain specific conditions (|97|129jl . The generalisation of 
section RHU to continuous mass, discrete time and arbitrary amount of mass transferred, 
yields a condition for factorisation with an appealingly simple general form. The 
factorisation property allows an exact analysis of steady state behaviour, in particular 
condensation. It should be noted that condensation is not restricted to models with 
factorised steady states, rather, factorisation, has offered, so far, the only opportunity 
to study condensation exactly. We remark that factorisation has also been used as a 
mean-field-type approximation to study condensation in systems where the steady state 
is not known exactly. 

More generally, it would be of interest to study the structure of non-factorised 
steady states and how condensation arises. For example, one generalisation of a 
factorised steady state is to a matrix product state. 

Returning to the case of factorised steady state, there have been a number of 
developments in our understanding of condensation. Here we have reviewed the distinct 
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mechanisms of condensation in homogeneous systems f sections EH I6.3.7|) . heterogeneous 
systems including the case of a single defect site 15.21 non-conserving systems 16.41 
and systems with more than one species of particle 16.11 We have also seen an example 
where the condensate is mobile rather than fixed 16.3.71 In particular the analysis of 
condensation within the canonical ensemble 14.21 should allow a deeper understanding of 
the nature of condensates. 

We have also reviewed the dynamics of condensation. The coarsening dynamics 
is now typically understood at a mean-field level or using heuristic random walk 
arguments. These arguments have been applied to condensation in homogeneous 14.51 
heterogeneous 15.31 and two-species systems 16.1.31 An open question is whether one can 
obtain scaling functions for dynamics beyond the mean field level. This would enable 
one to address questions of universality. We also note that for certain parameters (for 
which condensation does not occur) the Bethe ansatz has been used to calculate exactly 
dynamical properties of the ZRP [99 . 

Finally, we are aware that this review is of a theoretical nature and we have not 
emphasized or speculated upon experimental realisations. In section 13.11 we discussed 
clustering in shaken granular gases, which furnishes a pleasing experimental example 
of condensation in interacting particle systems. Another realisation of the ZRP noted 
in |lUUj concerns the exchange of monomers between protein filaments. It remains an 
important challenge to establish further experimental instances of interacting particle 
systems and the associated phenomena reviewed here. 
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